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CHAPTER  1 
INTRODUCTION 
For  nearly  two  decades,  there  has  been  a  large  and  continuing  interest 
in  the  study  of  maintenance  models  for  items  with  stochastic  failure.  This 
interest  has  its  roots  in  many  military  and  industrial  applications.  Lately, 
however,  new  applications  have  arisen  in  such  areas  as  health,  ecology,  and 
the  environment.  Although  it  is  not  possible  to  detail  these  many  appli- 
cations of  maintenance  models,  some  of  them  are:  the  maintenance  of  complex 
electronic  and/or  mechanical  equipment,  maintenance  of  human  body,  inspection 
and  control  of  pollutants  in  the  environment,  and  maintenance  of  ecological 
balance  in  populations  of  plants  and  animals.  This  interest  has  been  provoked 
by  the  high  cost  and  extraordinary  demands  made  of  modern  equipment  like  jet 
liners,  electronic  computers,  ballistic  missiles,  etc.  Operational  require- 
ments can  be  achieved  only  by  observing  relatively  sophisticated  maintenance 

The  practical  need  for  subtle  and  delicate  maintenance  policies  has 
stimulated  theoretical  interest  and  in  many  cases  has  led  to  the  development 
of  policies  that  posses  theoretical  novelty  and  practical  importance.  The 
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justified  when  the  cost  of  unscheduled  maintenance  action  is  higher  than 

the  cost  of  scheduled  maintenance.  The  name  of  preparedness  has  been  giv 
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to  its  policy  since  in  some  cases  there  are  standby  emergency  systems  which 
are  put  into  operation  when  the  main  system  suddenly  fails.  These  standby 
systems  should  be  checked  for  their  preparedness;  but  not  all  of  the  appli- 
cations of  preparedness  maintenance  modelling  lies  within  emergency  systems. 
There  has  been  already  a  lot  of  work  done  in  both  preventive  maintenance 
[1,2,4,5,10,12,19,23,28,31]  and  preparedness  maintenance  [1,3,8,11,14,18, 
23,28,29,32]  areas.  The  basic  criterion  is  to  minimize  the  total  cost  of 
maintenance  by  optimum  scheduling  of  the  maintenance  actions  such  as  repair, 
inspection  or  replacement. 

The  optimization  objective  for  preparedness  inspection  models  is  to 
optimally  balance  the  cost  of  undetected  system  failure  over  the  inspection 

Assuming: 

f(t)  »  The  density  function  of  the  time  to  failure  of  the 

I  »  The  cost  of  an  inspection  except  at  time  tQ  where  IQ=0; 

a  =  The  cost  per  unit  time  associated  with  an  undetected 
failed  system; 
and  that  the  inspections  to  be  performed  at  times  tj ,  t2>  t,,  ...,  until 
the  system  failure  is  detected  (See  Fig  1.),  then,  if  failure  occurs  between 
time  tQ  and  t, ,  say  at  X, ,  the  cost  of  the  cycle  of  operation  (See  Fig.  1.) 
would  be: 

KD  ♦  a(t,  -  x,)  (1.1) 

and  the  expected  value  of  this  cost  in  the  interval  [tQ,  t,]  is: 

/  [1(1)  +  a(t]  -  x)]  f(x)dx  .  (1.2) 

t0 

If  failure  occurs  between  t.  and  t9,  say  at  time  x,,  the  cost  of  the 
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Fig.  1.  The  Illustration  of  One  Time  Span  or  One  Cycl 


cycle  in  the  interval  [t^  t~]  would  be: 

1(1+1)  +  a(t2  -  x2)  (1.3) 

and  the  expected  value  of  this  cost  would  be: 

h 

J     [1(1+1)  +  a(t2  -  x)]  f(x)dx  .  (1.4) 

h 

In  a  manner  similar  to  the  above,  the  costs  and  probabilities  of  all 
possible  cycles  can  be  determined  to  give  the  expected  cost  per  cycle  as: 
Expected  cost  per  cycle 

=  /  [1(0+1)  +  a(t,  -  x)]  f(x)dx 


-  /  [1(1+1)  +  a(t2  -  x)]  f(x)dx 

t3 

-  /  [1(2+1)  +  a(t3  -  x)]  f(x)dx 


Therefore  if 

L(t-|,   t~>   t,,    . ..)  «  Expected  Cost  Per  Cycle 

then  Eq.    (1.5)   can  be  written  as 

k-n-1     Vl 
L(tr  t2,  t3.   ...)  -     I         /         [I(k+1)  +  a(tR+1   -  x)]  f(x)dx  . 

k=0    h  (i.6) 

n  -  1,  2,  3,  ....  - 

The  objective  function  to  be  minimized  is  Lltj,  t-,  t,,  ...)  if  the 
failed  system  is  replaced,  renewed,  repaired  etc.  at  the  end  of  its  fixed 


maximum  life  time.     If  the  system  is  replaced  or  renewed  or  repaired  when 
failed,  then,  the  objective  function  to  be  minimized  would  be  expected 
cost  per  unit  time  which  is: 

Expected  Cost  Per  Unit  Tire  .  Expected  Cost^Cycle   .  (w) 

The  optimum  solution  gives  the  optimum  number  of  inspections  n,  and 
optimum  timing  of  inspections,  i.e,   t, ,  t„,   t.,    ...,  t   .     Most  of  the 
literature  [6,   17,  20]  is  about  the  optimum  solution  which  minimizes  the 
total   maintenance  cost  per  unit  time  when  the  failure  distribution  of  the 
system  is  completely  known  in  advance.     But  there  are  many  cases  where  the 
failure  distribution  of  the  system  is  completely  unknown  or  is  partially 
known.      In  these  cases  different  methods  have  been  devised  depending  on 
the  kind  and  amount  of  information  which  is  available  about  the  failure 
characteristics  of  the  system  [26]. 

Minimax  policies  have  been  devised  to  cope  with  the  situation  in  which 
the  decision  maker  has  virtually  no  information  about  the  failure  distri- 
bution.    This  method  minimizes  the  maximum  possible  loss  or  maintenance  cost 
that  can  occur  due  to  the  failure  characteristics  by  designating  the  optimum 
number  of  inspection  and  optimum  timing  of  them.     If  P  =  {t,,  t?,  t,,   ...,  t  } 
is  the  space  of  the  inspection  times  t, ,   t«,  t,,    ...,   t  ,  and  if  F  is  any 
constraints  of  the  problem  related  to  the 
system,  then,   the  minimax  policy  can  be  formu- 

ixF  L(P.F),  (1.8) 

intenance  cost  per  cycle  as  given  by  Eq.  (1.6). 
expected  maintenance  cost,  the  optimum  number 
and  the  failure  distribution  which  gives  the 
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maximum  possible  rise  to  the  total  maintenance  cost  L  respectively.  The 
minimax  policy  is  not  limited  to  preparedness  or  inspection  models  and  has 
been  implemented  in  the  area  of  preventive  maintenance  by  Barlow  [4]. 
But  the  major  work  in  minimax  policy  has  been  done  for  the  preparedness  and 
inspection  policies.  In  minimax  policies  applied  to  inspection  models  the 
time  horizon  T  is  assumed  to  be  finite;  i.e.,  the  cost  accounting  stops  at 
either  the  first  inspection  to  detect  failure  or  at  time  T,    whichever 
happens  first.  The  reason  for  this  is  that  for  any  possible  inspection 
schedule  there  exists  a  distribution  which  would  induce  an  arbitrarily  high 
expected  cost  during  an  infinite  time  horizon.  Hence,  a  minimax  solution 
would  not  then  exist.  The  finite  horizon  assumption  can  be  found  to  have 
many  implications.  As  an  example,  consider  the  problem  of  detecting  the 
occurrence  of  an  event  (say,  the  arrival  of  an  enemy  missile  or  the  presence 
of  some  grave  illness  such  as  cancer)  when  the  time  of  occurrence  is  not 
known  in  advance.  Each  inspection  involves  a  cost  so  that  we  do  not  wish 
to  check  too  often.  On  the  other  hand,  there  is  a  penalty  cost  associated 
with  the  lapsed  time  between  occurrence  and  its  detection  so  that  we  wish 
to  check  often  enough  to  avoid  a  long  lapse  of  time  between  failure  and  its 

Derman  [9]  has  found  analytically  the  optimum  number  of  inspections  and 
timing,  applying  the  minimax  policy,  when  the  failure  distribution  is  completely 
unknown,  to  the  loss  function,  given  by  Eq.  (1.6),  under  the  following  assump- 

1.  The  maximum  life  or  service  time  of  the  system  T  is  limited  and 

known,  i.e.,  until  time  T  the  system  has  either  failed  or  the  system 


will  be  put  out  of  service  or  repaired,  renewed  etc.  at  the  end 
of  fixed  periods  of  length  T  equal  to  the  maximum  life  time  of 
the  system. 

2.  The  failure  can  only  be  detected  by  inspection  with  a  certain 
probability  P'(P'  >  0)  and  the  inspection  does  not  affect  the 
failure  characteristics. 

3.  The  inspection  time  is  negligible. 

4.  Each  check  entails  a  cost  I. 

5.  The  time  elapsed  between  system  failure  and  its  discovery  at  the 
next  check  has  a  cost  per  unit  of  time  a. 

In  many  cases  there  are  estimable  costs  associated  with  the  resumption 
in  service   (storage)  of  a  unit  which  has  failed.     If  the  unit  is  a  production 
system  the  costs  are  associated  with  the  amount  of  defective  product  produced, 


n  kind  of  missile  fuels)  the  costs 
implications  of  using,   unknowingly, 


if  it  is  material  in  storage  (e.g.,  cert 
are  derived  from  considering  the  various 
the  unserviceable  material;    ...   and  the  like. 

Derman  [9],  proved  that  under  the  above  assumptions  the  minimax  schedule 
is  given   by 

t     =  ,p.    f_E_+  j_(n[(n+l)P'+2l  ,   ,.      .]  fl     , 

Li  jnP'+l       2a   l       nP'+l  u    |;J  w's; 

1  =0,   1,    ...,   n, 

where  n,  the  number  of  inspections,  is  the  largest  integer  such  that 

IP,2n2  +  IP'(2-P')  n+2(I-P'aT)  <  0.  (1.10) 

The  minimax  expected  cost  L**,  when  P'=l  is  given  by 

,**  -  aT  ,1  n(n+3)  „  „, 


and  t.   is  the  time  of  ith  inspection.     Roeloffs   (30],  obtained  analytically 
the  minimax  schedule  for  the  German's  Problem  [9]  with  the  further  assumption 
that  the  location  x',  of  the  100. P      percentile  of  the  otherwise  unknown 
cumulative  failure  distribution  function  F,  of  the  system  is  known.     That  is, 

F(t-X')  =  p  X'    >  0;       0  <  p  <  1.  (1.12) 


Oit1i...itJ<xl<V]1...<Vn<T.  (1.13) 

The  minimax  inspection  schedule  and  the  minimax  expected  cost  analytically 
obtained  by  Roeloffs  [30]  does  not  have  a  simple  form  as  Derman's  does  and 
is  more  sophisticated  and  involved.     The  Roeloffs'    solution  is  also  under  the 
assumption  that  the  probability  of  the  detection  of  failure  p',  upon  in- 
spection is  one.     Roeloffs  also  showed  that  the  added  information  about  the 
location  of  a  percentile  of  the  failure  distribution  improves  the  minimax 
expected  loss  or  maintenance  cost  function  given  by  Eq.    (1.6)  and  results  in 
less  maintenance  cost  compared  to  the  solution  given  by  Derman,   i.e.,  Eqs. 
(1.9-11). 

The  implication  of  dynamic  programming  in  the  area  of  maintenance  is 
well   known  to  the  researchers  in  this  field.     Hasting  [13],  Jardine  [16,   17], 
Bellman  [7]  and  many  others  have  applied  the  dynamic  programming  to  a  variety 
of  replacement- repair  maintenance  problems,  but  to  the  best  of  my  knowledge, 
all  of  them  except  Kander  [20,21],   have  assumed  a  known  system  failure  dis- 
tribution and  in  some  cases  even  more  specifically  they  have  assumed  a 
certain  type  of  failure  distribution.     Kander  solved  several   different  kinds 


of  inspection  scheduling  problems  by  converting  the  minimization  of  the  loss 
function,  i.e.,  Eq.  (1.6),  into  optimization  of  recurrence  relationships  [21]. 

Introduction  of  the  implication  of  dynamic  programming  in  minimax  policy 
for  the  solution  of  the  problem  of  optimum  inspection  frequency  and  timing 
has  been  accomplished  by  Kander  [22],  however,  the  idea  and  theory  of  minimax 
dynamic  programming  was  well  introduced  and  established  by  Bellman  [7],  the 
pioneer  in  this  field  and  later  repeated  by  Jacobs  [15].  Kander  implemented 
the  dynamic  programming  methodology  and  combined  it  with  minimax  policy  to 
obtain  numerically  the  optimum  solution  to  the  Roeloffs'  problem,  discussed 
earlier,  with  a  further  assumption  that  the  system  failure  distribution  is 
IFR  (Increasing  Failure  Rate).  The  detail  of  this  method  has  been  fully 
explained  in  Chapter  2  since  the  present  work  is  highly  based  on  the  Kander's 
solution  procedure;  but  for  the  time  being  it  should  be  mentioned  that  Kander 
showed  that  his  solution  reduces  the  total  maintenance  cost  compared  to 
Roeloffs'  solution  [30]. 

The  purpose  of  the  present  work  has  been  directed  toward  obtaining  an 
upper  bound  for  the  expected  total  cost  of  the  system  maintenance  when  op- 
timum policy  is  applied  to  the  actual  failure  distribution  of  the  system 
and- under  Roeloffs'  assumptions  [30]  with  addition  to: 

A.  The  locations  x^  and  x£,  of  the  100. p1  th  and  100. p-  th  percentiles 
respectively  of  the  otherwise  unknown  life  or  failure  distribution 
of  the  system  are  known.  That  is 

F(t]  =  x})  »  p]      x^  >  0;     0  <  p1  <  1.   (1.14) 
F(t'  ■  x')  =  p2      x'  >  0;     0  <  p2  <  1.   (1.15) 


In  order  that  an  IFR  distribution  pass  through  these  two  points  it  i 
necessary  for  p1 ,  p2,  xj  and  x£  to  satisfy  the  following  inequality 


B.  The  inspection  times  can  only  occur  at  discrete  points  in  time 
between  time  0  to  T. 

Kander  did  not  consider  the  restriction  on  the  inspection  time  imposed 
by  assumption  8.  Assumption  B,  which  reduces  the  computation  time  when 
utilized  in  the  dynamic  program,  is  a  realistic  assumption  in  many  cases,  i.e 
in  cases  where  the  inspector  is  available  for  inspection  only  at  certain  time 
and  not  all  the  times.  Assumption  A,  which  is  different  from  Kander' s  by 
the  information  about  the  location  of  two  percentiles  of  the  system  failure 
distribution  instead  of  one,  is  utilized  to  give  the  relative  value  of  the 
added  information  and  the  reduction  in  estimation  of  the  total  inspection  or 
maintenance  cost  of  the  system. 

In  Chapter  2,  the  properties  of  IFR  distributions  in  which  the  basic 
formulation  lies  are  presented.  The  loss  function  or  Eq.  (1.6),  has  been 
rearranged  in  an  order  to  be  suitable  for  dynamic  programming  formulation. 
The  recurrence  relationships,  state  and  control  variables,  stages  and  logic 
of  the  problem  have  been  formed  and  explained  in  detail. 

In  Chapter  3,  the  details  of  the  computational  procedure  for  the 
computer  programming  of  the  two  models,  i.e.,  with  information  about  one 
point  and  two  points  of  the  system  failure  distribution  are  stated. 

In  Chapter  4,  the  convergence  and  accuracy  of  the  solutions  are  being 
shown  through  an  example  and  then  the  results  are  presented  in  the  form  of 


tables  and  figures  showing  the  sensitivity  of  the  upper  bound  for  optimum 
expected  total  system  maintenance  cost  !_**,  respect  to  one  and  two  known 
locations  of  the  percentiles  of  system  failure  distribution  together  with 
sensitivity  of  L**  respect  to  the  values  of  the  elements  of  the  information 
parameter  vector  which  consists  of  elements  like  a_,  the  cost  of  undetected 
failure,  and  l_,   the  cost  of  every  inspection.  The  value  of  information  v.s. 
the  location  of  the  known  percentiles  of  a  I FR  distribution  has  been  illus- 
trated by  an  example.  Also  an  example  about  the  application  of  the  models 
is  given. 

Chapter  5,  gives  the  conclusion  derived  from  the  outputs  of  the  compute 
program  and  results  of  Chapter  4  plus  the  comparison  of  the  upper  bound 
expected  total  system  maintenance  cost,  having  information  about  one  point 
of  the  failure  distribution  (Model  A)  and  two  point  of  the  failure  distri- 
bution (Model  B).  This  will  be  followed  by  the  possibilities  for  further 
research  in  this  area. 
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L*=L(PQ,  F0)  =MinL(P,F0)lL(P,F0),  (2.1.!) 

where  P  is  any  inspection  policy. 

The  method  devised  for  partial  knowledge  (as  detailed  above)  assumes 
a  distribution  F**  while  subsequently  deriving  an  optimal  policy  P**  with 
minimal  loss  L** 

L**  =  L(P**,F**)  =  Min  L(P,F**).  (2.1.2) 

P 

In  general  L*  can  be  larger  or  smaller  than  L**,  but  the  virtue  of 

minimax  policy  is  that  L**  represents  generally  an  upper  bound. 

For  a  minimax  policy  we  have 

L(P,F**)  »  Max  L(P,F)  >  L(P,F),  (2.1.3) 

F 

and      L**  =  L(P**,F**)  =  Min  L(P,F**)  =  Min  Max  L(P,F).        (2.1.4) 
P  P   F 

If  a  saddle  point  exists  then: 

L**  =  Max  L(P**,F). 

F 

Now  having  L**  defined  by  Eq.  (2.1.4)  we  can  write 

L**  •  Min  Max  L(P,F)  >L(P0,FQ)  ■  L*  (2.1.5) 

The  proof  is  as  follows. 
Proof:  If  (2.1.5)  is  not  true  then  we  should  have  L**  <  L*.  But  according 
to  (2.1.1)  we  have  L*  <  L(P,FQ)  for  any  policy  P  including  minimax  policy 
P**  so  that 


On  the  other  hand  according  to  (2.1.3)  we  have  L(P,F**)  >  L(P,F)  for 
any  policy  P  and  failure  distribution  F  including  P**  and  FQ  respectively 
so  that  we  can  write  L(P**,F**)  >  L(P**,FQ).  But  L(P**,F**)  according  to 
2.1.4  is  the  minimax  cost  L**  so  we  should  have  L**  >  L(P**,F  ).  This 
contradicts  Eq.  (2.1.6)  derived  on  the  assumption  that  L**  <  L*.  So  it  is 
proved  that  indeed  L**  >  L*.  That  is  the  minimax  loss  or  cost  L**  is  an 
upper  bound  for  the  optimum  inspection  policies  P  applied  to  failure  dis- 
tributions  FQ. 

In  this  work  improved  upper  bounds  for  the  optimum  total  expected  cost 
will  be  obtained  by  assuming  the  knowledge  about  the  location  of  one  and  two 
percentiles  of  f.d.  and  also  that  failure  distributions  are  IFR  (Increasing 
Failure  Rate),  using  dynamic  programming  methodology. 

In  this  chapter,  first  the  assumptions  of  the  present  study  are  stated, 
followed  by  the  mathematical  preparation,  model  presentation,  formulation  of 
the  loss  function  and  finally  the  dynamic  programming  formulation  of  the 

2.2.  THE  PROBLEM  STATEMENT 

The  objective  is  to  find  the  upper  bound  cost  for  the  optimum  inspection 
policies  for  known  failure  distributions  L**  under  the  following  assumptions: 

(a)  The  probability  of  the  failure  of  the  system  at  time  T  is  one. 
This  means  that  T  is  the  maximum  life  time  of  the  system  and  the 
inspection  ends  some  time  between  first  inspection  at  time  0  and 
time  T  upon  the  detection  of  failure. 

(b)  The  period  for  renewal,  replacement,  replenishment  etc.  is  fixed 
and  is  equal  to  T. 


(c)  A  system  failure  is  detected  only  through  inspection,  which  costs 
X  dollars  each.  Inspection  is  carried  also  at  times  t  ■  0  and 

t  =  T. 

(d)  The  time  elapsed  between  system  failure  and  its  detection  at  the 
next  inspection  costs  a_  dollars  per  unit  time. 

(e)  Inspection  takes  negligible  time,  the  system  cannot  fail  during 
an  inspection  and  is  not  degraded  by  inspection. 

(f)  The  inspection  can  only  be  performed  at  certain  discrete  points 
in  time. 

(g)  The  failure  distribution  of  the  system  is  IFR  (Increasing  Failure 
Rate). 

The  minimax  solution  gives  the  upper  bound  for  optimum  expected  total 
maintenance  cost  per  cycle  L**  according  to  (2.1.5)  where  L  is  given  by 
Expected  Maintenance  Cost  Per  Cycle,  L 

n-1         rn-1  *„  , 

-   I  I  F*^)  «■  a  \\     F*(t  )(t    -  t  )  -  /   F*(t)  dt  ,    (2.2.1) 
i-0         1.1*0    '   1  '    1   0        J 

and  n  is  the  total  number  of  inspections  and  F*(t)  ■  1  -  F(t)  is  the 
probability  that  the  system  has  not  failed  until  time  t. 

The  failure  distribution  F(t)  is  unknown  except  at  the  time  T  where 
F(T)  ■  1  and  also  one  of  the  following  two  additional  informations  depending 
on  one  of  the  two  models  A  or  B  is  given.  That  is 

For  Model  A 

The  location  x',  of  the  100. P  th  percentile  of  the  otherwise  unknown 
failure  distribution  of  the  system  is  known.  That  is 

F(t'  -  x1)  =  p      x1  >  0;     0  <  p  <  1.  (2.2.2) 


For  Model   B 

The  locations  xj  and  x£,  of  the  100. p]   th  and  100. p2  th  percentiles  re- 
spectively of  the  otherwise  unknown  failure  distribution  of  the  system  are 
known.     That  is: 

F(t'    -  xp-  p1  x\  >  0;  0  <  Pl   <  1,  (2.2.3) 

and  f(f2  =  xp  =  p2  x2>0;  0  <  p2  <  1 .  (2.2.4) 

In  order  that  IFR  distributions  pass  through  these  two  points,   i.e., 
(xj,  p-,)  and  (x2,  p2),   it  is  necessary  for  p] ,  p2>   xJ   and  x2  to  satisfy  the 
following  inequality 

1-Po 

Log   (1-Pl)  Log  (j^) 

"— *? -"      x'    -   x'1  <2-2-5' 


2.3  MATHEMATICAL  PREPARATIONS 

According  to  [6],  a  failure  distribution  F(u),  u  >  0  which  is  IFR, 
crosses  the  exponential  distribution  l-e"aU,  a  >   0  from  below  at  most  once 
in  addition  to  coincidence  at  u  =  0  and  =°  unless  they  coincide  identically. 
This  can  be  easily  seen,  that  is  if  F  is  any  IFR  and  the  exponential  dis- 
tribution is  l-e"aU  then  the  two  functions  agree  at  the  origin  (See  Fig.  2) 
and  equating  the  two  functions 

F(t)  =  l-e"at  (2.3.1) 

and  shifting  F(t)  to  the  right  and  e"at  to  the  left  of  the  equal  sign  and 
calling  l-F(t)  as  F(t)  and  extracting  the  logarithms  then  : 

G(t)  ■  log  F(t)  -  (-at).  (2.3.2) 


Fig.  2.  Exponential  and  Arbitrary  IFR  Distribution,  F(t) 
Crossing  Exponential  Distribution  from  Below 


It  is  clear  that  G(t)   is  zero  at  t  =  0  since  for  t  =  0,  F(t)  =  1  and 
logel  =  0,  so  G(t)  =0-0=0.     According  to  the  definition  of  failure  rate 

Failure  Rate  =  r(u)  =     F<u+&)   '  F  <u>  (2.3.3) 

F(u) 

where  a  >  0  is  an  increment  of  time.     The  definition  of  an  IFR  distribution 

signifies  that 

p(u2)>r(Ul),         u2>ur  (2.3.4) 

Applying  Eq.    (2.3.3)   to  exponential    distribution  1   -  e"at, 

,      -a(t+A)         M      -at. 

Exponential    F.R.   =  r  (t)  =  ■L-^ K'~e L  =  i   .  e_aA 

l-(l-e-at) 

=  constant,   (2.3.5) 

and  calling  the  failure  rate  of  the  arbitrary  IFR  distribution  F(t)  as  r(t), 
then  F(t)  crosses  1  -  e"at  at  most  one  more  time  from  below  (See  Fig.  2.). 
F(t)  cannot  cross  1  -  e"a  from  above  since  if  that  happens  then 
r(t*)  <  re(t*)  =  1  -  e'aA  =   constant  where  tf  >  0  is  the  first  crossing  point 
after  the  origin.  But  the  first  crossing  from  above  cannot  happen  unless 
r(0)  >  re(0)  =  1  -  e'^constant.  This  means  that  r(0)  >  1  -  e"aA  >  r(t*) 
which  violates  the  condition  for  F(t)  to  be  an  IFR.  As  a  result  the  first 
crossing  of  1  -  e    by  F(t)  happens  from  below.  The  proof  that  F(t)  cannot 
cross  1  -  e"011-  more  than  once  is  similar  to  the  previous  one,  i.e.,  if  F(t) 
crosses  1  -  e    for  the  second  time,  it  should  be  from  above  which  means 
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first  and  second  crossing  points  besides  origin.  But  F(t)  has  already  crossed 
the  exponential  distribution  from  below  (proved  previously)  at  t?  so  we  already 
have  r(t*)  >  rg(t*)  =  1  -  e"aA=  const.  This  means  that  r(tf)  >  1  -  e"aA 
>  r(tJi)  which  again  violates  the  condition  for  F(t)  to  be  an  IFR.  So  no  IFR 
distribution  can  cross  exponential  distribution  more  than  once  besides  at 
t  =  0  and  t  =  »  and  this  crossing  should  take  place  from  below.  Consequently 
for  F(u  =  x)  =  p  the  common  point,  we  obtain: 


F(u)  <  1  ■ 


if  u  < 


(2.3.6) 
(2.3.7) 


introduce  a  transformed  exponential  distribution  function  - 


from  Eq.  (2.3.3)  it  follows  that  the  transformed  distribution  is  of  constant 


r(v)  =  1  -  e  ",  u  >  a,   A  >  0  .  (2.3.10) 

The  transformed  exponential  distribution  has  the  following  properties  [22]: 


PROPERTY  1:  Generalizing  properties  of  the  exponential  distribution 
(e.d.),  the  t.e.d.  can  be  crossed  by  an  IFR  distribution  at  most  twice,  the 
first  taking  place  from  above. 

REASON:  Assume  that  F(t),  which  is  IFR,  crosses  the  t.e.d.  F(v),  with 
constant  failure  rate  r(v),  first  from  above  and  then  from  below.  At  the 
first  coincidence  r(t1)  <  r(v)  and  at  the  second  r(t2)  >  r(v)  so  that  F(t) 
can  indeed  be  IFR.  A  further  crossing  from  above  would  mean  r(t  )  <  r(v), 
contradicting  the  assumption  of  IFR  distribution. 

Further,  we  can  uniquely  find  parameters  c,  d  such  that  the  t.e.d. 
passes  through  given  two  points  (t, ,  p. ) :  p.  =  F(t  =  t. ) ,  i  =  1 ,  2. 

PROPERTY  2:  Given  three  points 

(tr  Pl):  p.  =  F(w  =  ti),   1-0,1,2   tQ  <  t,  <  t2.  (2.3.11) 

(a)  An  IFR  distribution  F(w)  passes  through  the  three  points  only  if 
(tr  Pl)  lies  on  or  below  the  t.e.d.  through  (tQ,  pQ)  ,  (t2,  p2). 

(b)  The  IFR  distribution  F(w)  which  passes  through  the  three  points  and 

h 

possesses  maximum  area  /   F(w)  d  ,  is  given  by  the  two  t.e.d. 's  which  meet 


-  e  '  (2.3.12) 

-TTT^  Loge  (-rr^1)  (2-3.13) 

-  Loge  (1  -  p.)  +  t.  (2.3.14) 

w  it,   i  =  1,  2,  (2.3.15) 


c2  >  c 


REASON  (a):  This  property  is  a  direct  result  of  property  1. 
REASON  (b):  Let  us  define  the  following  distribution  functions  (See 
Fig.  3.):  t.e.d.  connecting  points  (t2>  p2)  ,  (tQ,  pQ) 

F3(w)  -  1  -  e" 

t.e.d.   connecting  points   (t,,   P-,),    k*.q,   Vq, 

w  -  d  ; 
F](w)  =  1   -  e     '  ' 

t.e.d.  connecting  points  (t2>  p2), 

F2(w)  «  1  -  e  °' 
We  assume  according  to  (a)  that 


Since  F^w)  crosses  t.e.d.  F3(w)  from  above  (See  Fig.  3.)  it  folic 
jctly  from  definition  (2.3.3)  that: 


r3(t0)>r1(t0),  (2.3.21) 

where  r.(-)  is  the  failure  rate  for  Fi ( • ) .  Similarly  since  F2(w)  crosses 
F3(w)  from  below  (See  Fig.  3.) 

r3(t2)  <  r2(t2),  (2.3.22) 
so  that  from  (2.3.21),  (2.3.22)  we  obtain: 

^(ViW  (2"3-23) 


from  which  by  (2.3.10)  we  have: 


The  Eqs.  (2.3.12  -  15)  can  easily  be  obtained  since  F^w)  and  F,(w)  have 
point  (t,,  p-,)  in  common. 

Viewing  now  F,  (w),  w  <_  t,  and  F„(w)  >_  t,  as  parts  of  one  distribution 
F(w),  we  obtain  for  the  failure  rate  of  the  latter  from  (2.3.3),  (2.3.10) 
and  (2.3.12  -  16): 


:2(w+A-  d2)       -c2(a-  d2)   ,  _c  . 


which  is  increasing  in  w,  so  that  F(w)  is  IFR. 

From  Property  1  it  follows  that  no  IFR  distribution  can  pass  through 
points  (t,,  p.),  i  =  1,  2  and  above  F(w).  Therefore  indeed 


/       F(w)  c 


n  among  all    IFR  dis 


2.4     THE  OBJECTIVE   FUNCTION 

In  Chapter  1    it  was  shown  that  the  ex 
cycle  could  be  found  from  Eq.    (1.6)  which 


n  =  1,  2,  3,  ....  ». 
This  equation  is  valid  when  there  is  no  maximum  life  time  for  the  system 
and  inspection  continues  until  the  system  fails.  According  to  the  assumption 
a  of  Section  2.2  there  is  a  maximum  life  time  T  for  the  system,  i.e., 
F(T)  ■  1  for  both  model  A  and  model  B  presented  here.  Also  the  definition  of 
the  problem  signifies  that  inspection  is  performed  at  time  tQ  =  0  and 
t  =  T.  With  these  points  in  mind  Eq.  (1.6)  can  be  written  as 


^=n-2 


L(t  =  0,  t   t2,  ...,  t  -T)  -  I         /    [1(1*1)  ♦  a(t.+,-x) 
k=0   tk  k+l 

t  "T 

]  f(x)  dx  +  /  "  [I(n)  *   a(T-x) 
tn-l 

]  f(x)  dx,  (2. 

re  N  is  the  maximum  possible  number  of  inspections  between  t  =  0  to 
T  including. 
Equation  (2.4.1)  can  be  written  in  the  following  form 

l  - 1  J    F*<V  +  a  fl  F*<V(Vi  -  h]  -  / "  F*ct)dt], 


(2.4.2) 


where  F*(t)  «  1  -  F(t)  and  F(t)  is  the  cumulative  failure  distribution  of  the 
system.  The  algebraic  manipulation  which  transforms  Eq.  (2.4.1)  into 
Eq.  (2.4.2)  is  presented  in  detail  in  Appendix  A. 


s  work  to  maximize  objective  fui 
on  policy  p  «  {tQ  ■  0,  t,,  tj, 

ons  F  which  pass  through  one  or 
i.e.,  A  or  B  respectively, 
butions  F(w)  leading  to  Max.  L 
sformed  exponential  distributioi 


2.5  FORMULATION 

It  is  the  objective  of  th 
by  (2.4.2)  for  a  given  inspect 
by  searching  all  IFR  distribut 
points  depending  on  the  model 

The  form  of  the  IFR  dist 
according  to  Property  2  by  tr 
(t.e.d)  connecting  all  points  (t. 
known  points. 
2.5.1  Model  A 

In  this  model  we  assume  that  one  point  of  the  failure  distr 
(x,p)  is  known.  If  m  is  the  number  of  inspections  before  x,  we 
time  sequence  tp,  t^ ,  t^,  ...,  t^,  x,  I  . ,  ...,  t  ,  according  X 
above  points  are  arranged.  We  recode  this  sequence  as 


■  F(t.)), 


■  0,    1,   2, 


(i)'   "(■ 


(x,  P)  ■  (1 


(v  pn) : 

tfing  these  we  cai 


o)=(t(0),P(Q)) 

m+D'  "(iTH-l)5,  c  =  c(rn 

=   (t(n)'%)>         Vl   =c 


functions 
,  n  and  the 


,  _t-cu>lw-d(1)' 


L(i-1)  i»il(1 

'(n-1)  -"-^r 


Log  (1  -  p) 


"(1) 


(1)  1C(2)  1  •■•  ±c(m+1)  i  ■•■  <C{ 

L0go  (V^-l) 


U(1) 


U(1) 


u(i)  ~  L(i-D 

-Log  (1  -  p   )  a 


u(i) 


p(i: 


J;i; 


u(1)' 


We  observe  that  F(w)  passes  through  points  (0,0),  (x.p)  and  possesse; 
jump  from  (t^,  p^)  to  tt*., ,  1)  (as  allowed  for  an  IFR  distribution  ir 
[6,  22]). 

The  objective  function  L,  given  feasible  ooints  (t.,  p.),  i  =  1,  2; 
■-.,  n-1,  is  indeed  maximized  by  function  F(w)  since 


/   F*(w)  dw  becomes  rr 


/   F*(w)  dw  -  /  "  (1  -  f(w))  dw 


tn=T 
•  T  -  /    F(w)  c 


but  /    F(w)  dw  becomes  maximized  according  to  Property  2  (2.3)  if  F(w)  is 
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given  by  the  set  of  the  transformed  exponential  distributions  defined  in 
(2.5.2).  The  value  of  maximum  life  time  T  is  fixed  so  Eq.  (2.5.4)  gives  the 
minimum  value  of  the  expression  for  all  IFR  distributions  F(t). 

The  feasible  region  of  F(t)  v.s.  t  diagram  is  shown  in  Fig.  4.  Between 
t  =  0  and  t  »  x  the  region  is  the  area  confined  from  above  by  F(t)  =  1  -  e"at, 
from  below  by  F(t)  =  0  and  from  left  and  right  by  t  =  0  and  t  =  x  respectively. 
The  reason  for  this  is  given  by  Property  (2-a-2.3)  which  states  that  any  IFR 
distribution  will  pass  through  three  points  (e.g.  [t  ■  0,  0],  [t-,  p.], 
[t  =  x,  p])  only  if  (ti ,  p..)  lies  on  or  below  the  t.e.d  through  the  first  and 
third  points  (e.g.  [t  =  0,  0],  [t  =  x,  p]).  Between  t  =  x  and  t  =  T  the 
region  is  the  area  confined  from  above  by  F(t)  »  1,  from  below  by 
F(t)  =  1  -  e"a  and  from  left  and  right  by  t  =  x  and  t  =  T  respectively  (See 
Fig.  4.).  The  reason  that  {t-,   p.)  cannot  be  below  F(t)  -  1  -  e"at  is  that 
any  IFR  distribution  which  crosses  F(t)  at  point  (t  =  x,  p)  should  have  a 
higher  failure  rate  at  this  point  than  1  -  e"4a.  So  no  point  (t-,  p.)  for 
ti  >  x  can  lie  below  F(t)  =  1  -  e"at.  In  addition  to  this  the  points 
(t. ,  p.),   must  satisfy  the  equations  given  in  (2.5.2). 


\\\\ 

, 

V\\\ 

F(t)  -  1-e 

.< 

A\\\\v 

l(t-x,p) 

Loge(l- 

) 

\ 

The  Feasible  Region  for  Points  (t^p^  for  Model  / 


The  following  relationship  exists  between  two  successive  points  {£*,   p^). 
(t._},   p._,).  That  is 

F*uciy  -*(t(w)>  •"Ct1,(tt1) " w 


This  can  be  obtained  from  equations  (2.5.2).  That  is 


F(t    v.,   .W(1-1)  "  u(i); 

since  both  points  (t/^.i)>  p(i-i)'  and  1^]'  P(i))  are  tne  two  ends  of  F(t)- 
Now  shifting  1  to  the  left  side  of  both  equations  and  multiplying  both  sides 
by  -1  and  dividing  by  each  other  then  we  have 


"-(D^(i)  "  u(1)  " 


lL(i-U' 

Maximization  for  given  policy  can  now  be  carried  out  by  dynamic  pro- 
graming methodology.  The  search  extends  over  all  p.  =  F(t.),  i  »  1,  2,  ...,  n. 
We  define  J.,   as  loss  or  maintenance  cost  on  the  interval  between  (i-l)th  and 
(i)th  inspection  [t^(  tj],  which  is  also  a  function  of  the  failure  prob- 
ability  at  (i-l)th  and  i  th  inspection.  That  is 


Ji  =  Vl-V  Pi-r  V  Pi  >  =  Ci  +  a(t.  -  t.^)]  F*(t1_l) 

-  a  /         F*(w)  dw  (2.5.6) 

ti-l 

ti   <  w  <  ti_1,     1  =  1,  2,   ...,  n. 
If  J?  =  Max     J.,  we  obtain  from  equations   (2.5.2)   that 

JT =  Ji(ti-r  pi-r  V  Pi'  =  Max   Ji  (^.v  p.-.i.  V  P  )        (2.5.7) 


m,m+2,...,n-l 


[I  +  a(tm+1  -  t|n)]F*(tm)  4 
F*(tm+1)   -    (1-p). 


;i   +  a(tn  -   tn-1)]  F*(tn.1)  if     1-n 

The  first  expressions  on  the  right  hand  side  of  equations   in   (2.5.7), 
i.e.,   [I  +  a(ti   -  t.^)]  F*(ti_1)   are  obviously  the  share  of  J. 
from  the  first  expression  in  the  expected  total  maintenance  cost  L  given 

n-1  n-1 

by  Eq.    (2.4.2),   i.e.,   I     j>     F*(ti)  +  a    J     F*(ti)(ti+1   -  tj).     The  second 


expression  a  /    F*(w)  dw,  is  also  the  share  of  J.  from  L,  i.e.,  a  /    F*(w)dw. 

*1-1  t  ° 

For  values  of  i  =  1 ,2,  . . . ,  m,  m  +  2,  . . . ,  n-1 ,  /   F*(w)  dw  can  be  calculated 
from  equations  (2.5.2).  That  is  i~1 


.  ^^j-i)  -  F*<y 

ci 
For  i  =  m+1 ,  F(w)  passes  through  the  known  point  (x,p),  sc 

/   F*(w)  dw  =  /*  F*(w)  dw  +  /  m+  F*(w)  dw 
•  /  e         dw  +  / 


:m(tm  "  dJ         "cm(x  "   V 


m+r"       m+V 


but  e     m  m     =  »    m+1  m+1     «  F*(x)  =   (1-p)  since  all   describe  the 

same  point  (x,p).     Also  cm  =  c  according  to  recodification  in  Section  5  of 
this  chapter.     So 


tm  cm+l 

For  i  =  n,   i.e.,   to  find  J*  between  one  to  the  last  inspection  at  t 
and  last  inspection  at  tn  =  T,  according  to  Eqs.    (2.5.7),   F(w)  =  1   or 
equivalents  F*(w)  =  0  so 
t  =T 
/         F*(w)   dw  -  0.  (2.5. 

Let  us  define  ICj(t.j.  Pi ,  c]  +  ])   as  loss  or  expected  maintenance  cost 
during  time   (0,   tj]  with  F(t,)  «  pi  when  the  parameter  of  the  transformed 
failure  distribution  connecting  point  (t,,   pf)   to   (t.+r  p .+] )   is  ci+r 
Then  we  have 


K,(tr  Pr   c2)  -  i}[tQ,   pQ,   tr  P])  -J^O,  0,   tr   Pl) 

Ki(ti'  pi'   ci+l]  =  Ki-l(ti-T   Pi-v   ci>  +  Wl'  Pi-r   V  pi> 

a  j     Jj^l-V   Pi_v  tv  Pj).  (2.5.11) 


Let  K  (T)  signify  the  loss  for  the  entire  time  period  T  at  F(T)  = 
any  t.e.d.   parameter  c     , 


Kn(T)  =  Kp(tn  -T,  P||-  1,  cn+1) 

=  Kn-l(tn-T  Pn-T   cn)  +  Wl*   pn-l  *   V   pn> 
=  L.  (2.5.12) 

From  this   it  should  be  clear  that  there  are  three  state  variables  in 
the  system,   i.e.,   ti ,  pi   and  c.+..     Optimization  can  now  be  carried  out 
in  two  phases.      In  the  first  phase  ^(t^,  Pm+1 ,   c^),   i.e.,  the  loss 
between  t  ■  0  to  the  time  of  the  first  inspection  larger  or  equal    to  x, 
tm+1  will   be  optimized  for  all   possible  values  of  t     .   and  p  +. .     In  the 
second  and  final   phase  Kn(tn=m+j  -  T,   pn=m+j  -  1,  cn+1=mfj+1   ■  -)  =  L  will 
be  optimized.     In  the  second  phase,  the  optimization  will   be  achieved,  taking 
the  optimum  Km+]   values  for  state  variables  at  time  t  -  t(n+1   as  the  K  values 
for  the  first  stage  of  the  second  phase  and  the  values  of  L  will   be  calculated 
for  all   possible  number  of  stages  both  in  phase  1   and  phase  2  and  the  optimum 
value  of  L  is  obtained  by  search  among  these  values. 

First  optimization  phase  -  Let  the  minimax  optimized  loss  for  time 
period   (0,   tj]  at  F(tj)  =  p.  when  the  parameter  of  the  transformed  failure 
distribution  connecting  (t^   p.)  to   (ti+] ,   P1+])   is  ci+]    be 

K${tv   p.,   ci+1)  =  Min  Max  ^.(t^   p.,   z.)    ,  (2.5.13) 

where  P   is  the  inspection  policy  and  F  is  the  failure  distribution  given  by 
Eqs.    (2.5.2).     By  dynamic  programming  procedure  we  then  obtain 


Kg(tp-0,  P0=0,  c0+1)=0. 
K*(tr   pr   C2)  =  J*(0,   0,   t,,   p1)  - 


1    ■   1 


p"(t-, )  <  p1  <  p'Ct,) 

K*(t2,   p2,  c3)  =  Min  Max{K|(tr  p] ,   Cg)  +  J*(t1 ,   p] ,   tg,  P2)} 

0<t2<x  0<tl<t2  V(t2>p2,c3,t1)<p1<S(t2,P2,t1) 

p"(t2)<p2<p'(t2)  i   =  2 


P"(t,)  <  p,  <  p'(t.) 


•  Min  Max{K^](tl..r   p^,  c,)  +  JfCt^-, ,  p._r   t.,  p.)}, 

0<ti_1<ti     v(tl.,Pi,ci.  +  1,t._1)<p._1<G(ti.,pi,t._1) 


K£tl(W    Pm+T    <=„,,) 


'"'Vi'^i^i'Vi1 


Max{K*(t  ,p   .c)  +  J*     (t, 


W^VVWW' 


lVT"! 


■tUVW 


0<Pm<G(tm) 

m  ■  0,1,2,    ...,   N, 


N]  is  the  maximum  possible  number  of  inspections  in  the  interval  (0,x) 
and  we  have  c^c^  •••<c<cm+1  since  the  failure  distribution  F(w)  is  IFR. 
Also  since  at  point  (x,p),  F(w)  crosses  the  exponential  distribution  F(u) 
(2.3.6-8),  we  should  have  o<c<c|n+1 ,  as  stated  already  in  Eqs.  (2.5.2). 
p'(tj)  and  p"(t.,)  which  are  the  upper  and  lower  limits  for  the  possible 
values  of  p1  in  the  feasible  region  (See  Fig.  5.)  are  given  by 


The  Feasible  Region  for  IFR  to 
Pass  Through  Points  (0,0)  and 
(x.p) 

The  Feasible  Region  for  IFR  to 
Pass  Through  Points  (0,0), 
(t.  ,,p.  ,)  and  (t. ,p.)  with 


Fig.  5.  The  Upper  and  L 


p'(x)  =  p"(x)  =  p  (2.5.17) 

P"(T)  =  1.  (2.5.18) 

In  order  for  any  IFR  distribution  to  pass  through  points  (0,  0), 
1*1-1*  pi-l'  and  1*1*  pi''  U  is  necessary  f°r  Point  (t._-|,  p._-|)  to  be 
located  below  the  transformed  failure  distribution  (t.e.d.)  through  points 
(0,  0)  and  (ti;,  p^ )  according  to  Property  (2-a-2.3)  for  t.  in  the  interval 
(0,x)  (See  t.e.d.  No  2  in  Fig.  5.).  The  parameter  of  the  t.e.d.  No  2  called 
cU.j.p.j)  can  be  calculated  from 


and  the  maximum  possible  value  of  p.  ,  given  by  G(t.,p-,t.  ,)  a 
to  Eq.  (2.5.5)  is 


G(tm)  is  the  particular  value  of  G  function  at  t_  and  can  be  found  with 

a  similar  type  of  reasoning  as  for  Eqs.  (2.5.19-20)  to  be 

-at 
G(tm)  ■  1-e   m.  (2.5.21) 

In  order  for  a  transformed  exponential   distribution  to  connect  points 
( t^ _-| » P-j _ -j )  ancl  (t^  ,pi )  and  have  a  parameter  c..   so  that  c.<c.  +  ,,  where  c.    . 
is  the  parameter  of  the  transformed  exponential   distribution  connecting 
points   (t. ,pi )  and  (t.+1>pi+1)  as  required  for  an  IFR  distribution   (See 
Fig.   5  t.e.d.   No  1),   the  point  (ti_1  ,pi_1 )   should  be  above  t.e.d.   No  1.     The 
reason  is  obvious.     Let  us  call    the  t.e.d.    passing  through  points   (t.,p.)  and 


(t.+1,p.+1)    as  v.,(t)  then  if  point  (ti_1,p._1)  is  located  b 
should  have  VjUi  ^)  >  p.,,.     Now  we  have 


)arameter  of  v,  (t)  = 


:an  write        6,-c  =  f^—  (Log^^-^)  -  Log^.p^     })) 


i   loge(- 


1-F(ti_1) 


since  V-j C*^_-|)  »  P^_-j   according  to  the  assumption  and  so  1-F  (t,  ,)  <  1   -  p. 

which  shows  that  loge(A<l).  But  log  (A<1)  is  always  negative.  This  means 
that  ci  -  c  =  B  where  B  is  a  positive  amount  so  ci  >  c  =  ci+r  But  we  shoul 
have  ci  £  ci+]  in  order  that  the  failure  distribution  to  be  IFR.  As  a  resul 
point  (ti_1  .Pn-_-j )  should  be  above  v,(t).  The  lower  limit  on  p.  is  given 
by  v  function  which  its  value  is  determined  by  the  above  mentioned  v,(t)  dis 
tribution  or  regional  boundaries  and  has  different  forms  at  different  stages 
and  different  points.  v(^>p1,ci+l't1-^  accord1n9  to  Eq.  (2.5.5)  is  given 
(t.-t.    ,) 


vCt,,p,,c,+ 


)  •  Max{l.   -   (1-p^e  i  +  1      '      '"'     or  0.0}    . 

(2.5.25) 


^Vl^im-l'V  1s  the  Particular  value  of  v  function  at  t+],  i.e.,  the 
lue  at  (rrt-l)th  stage  and  is  given  by 

"(WmH'V  "  Max{l-D-P)e  ""     or  0.0},     (2.5.26) 


The  value  of  c^i.e.,  the  parameter  for  the  t.e.d.  connecting  poin 

(ti,pi)  and  (t.^.p.^)  is  given  according  to  Eq.  (2.5.5)  by 

1-P,- 
Loge(T^T) 


We  now  choose  for  each  (tm+1.Pm+1)  that  m=m*  which  renders  K*   a  min- 
mum.  That  is 

K'WWV+rW  "  M1n{WWWVi)J 

*iV+i<T 

P"(tm*+l)<fV+l<p'(V+l)  m  =  0,1,2,    ...,  N1 

(2.5.29) 

C^V+rV^V+l^'V+l'V*0) 

■  Min  {R*         (T,l,»)}  ■  -  0.1,2,   ....  N,. 

1         (2.5.30) 

<*=m+1(T,l,<»),   the  total   expected  loss  when  there  is  only  one  inspection 

i  the  interval    (x,T]  is  given  by  the  last  equation  in   (2.5.14).     K**  „       is 


the  minima*  expected  total  maintenance  cost  in  the  time  interval  (0,tm*+1) 
where  x<t  *+1<T.  m*  is  the  optimum  number  of  inspections  in  the  interval 
(0,x).  From  the  above  formulations  it  becomes  clear  that  the  control  vari- 
ables are  the  previous  inspection  time  tj  ^  and  previous  cumulative  failure 
probability  p.  , .  The  stage  number  m  represents  the  number  of  inspections. 
Second  optimization  phase  -  we  continue  the  iterative  procedure,  taking  the 
K**  ]t+.  values  at  each  point  (t m*+1>P m*+i )  as  the  minimax  expected  maintenance 
cost  up  to  that  point  or  state  for  the  first  value  of  K**  at  the  first  stage 
in  the  second  phase  of  the  optimization.  That  is 

K**m*+j (tm*+j 'pm*+j ,cm*+j+l } 

P"(V+j)<Pm*+j<P'(v+j) 

Max{K*Wl(Wl,pm,+j.rW 

xiV+j.llV+j      v(  v+j.P|n,+j,v+J+1 ,  v+j)<pm,+j.1< 

<G(V+j,pm,+J.,V+j.1) 
+  WWrPm^rWW}>  (2.5.31) 

j   -  2,3,4,    ...,   N- 


=  Min  ^^WlWrWr"5 


j  =  2,3,    ...,  N2 

(2.5.32) 
where  Ng  is  the  maximum  possible  number  of  inspections  or  equivalently 
maximum  number  of  stages  in  the  second  phase.     Also  we  should  have 
a  1  cm*+1  1  cm*+2  £  •••  1  cn_i   according  to  the  definition  of  the  failure 
distribution  F(w)  given  in  Eqs.    (2.5.2-3).     P'Ct^J  and  P"(tm,+.)  are 
given  in  Eqs.    (2.5.15-18).     Function  G  with  a  similar  reasoning  as  for  phase 
one  is  given  by 

X) 


L°9j-lV 


n*+j^m*+j;  ~  ~       t 


Function  v  with  a  similar  reasonir 
Eqs.  (2.5.22-24)  is  given  by 

v(WWV+j+rW 

*  Max{l   -   (l-pm*+.)e  m*+j+1 

(2.5.34) 
and  function  v(t   .  J  is  a  particular  value  of  v  function  att.  .,. 


Wi*-1 


m*+j-l 


The  value  of  the  parameter  of  the  transformed  exponential  distributio 
which  passes  through  (tm*+J._-, .  Pm*+J-_-] )  and  It^y  Pm*+j)>  I.e.,  c^.  is 
given  by 


Log(- 


•m*+J 


K**n=m*+j(T,l,-),   given  by  Eq.    (2.5.32)   is  the  minimax  total   expected 
cost  in  the  interval   [0,T]  when  m*  inspections  are  performed  in  the  in- 
terval   [0,x)  and  j   inspection  performed  in  the  interval    [x,T].     Now  the 


hJ(T,l,-)  or 

j  ■  2,3,4, 


,(T, !,»)},  (2.5,37) 


where  £**n=rr]*+-|  (T.l  ■")   is  the  minimax  total   expected  cost  in  the  interval 
[0,T]  with  m*  inspections  in  the  interval    [0,x)  and  one  inspection  at  tims 
T  and  is  given  by  Eq.    (2.5.30).     j*  is  the  optimal   number  of  inspections 
in  the  interval   [x,T]. 

The  policy  p**  ■  {t|,  t|,  ...,  t*}  and  the  failure  distribution  F** 
defined  by  (2.1.3-4)  can  also  be  found  if  needed  by  a  backward  recursive 
procedure  from  the  computer  printout  (as  explained  in  Chapter  3). 


2.5.2  Model  B 

In  this  model  we  assume  that  the  location  x]  and  x2  of  the  lOO.p-,  th 
and  100. p„  th  percentiles  respectively  of  the  otherwise  unknown  failure  dis- 
tribution of  the  system  are  known.  If  m1  is  the  number  of  checks  or  in- 


spections before 
the  time  sequence 


;  the  number  of  inspectic 


according  to  which  above  points 
{(tm,  p.  .);  i  =  C 


-anged.  We  recode  this  seque 


t0,  pQ)  .  (0,  0)  =  (t(0), 

W 

-„'-!'   Pn'.-|)  =   (t(n')'   p 

n')>       c 

'  =  C(j^+1) 

<2'   h]  =   ^(VZ)'   p(m2+- 

,).         = 

=  °(m2+2) 

-n-l-  Pn_!)  "   tt(n+1).   P( 

♦1>>       c 

-1   =  c(n+l) 

=  obtair 


i  define  F(w)  with  a  series  of  transformed  exponent! 


l-oge(  l-p1 ) 


-Pi  _     LogeO-Pi) 


C(D  -  •••  ^(m,)  -  c(|^+l)  -  •■•  i  c(m2+1)  -  c(m2+2)  ±  —  ±  V 
al   "  C(m1+1)   *  a2  c(m2+2)  *  a2 


(D  ~  tM,  -  tn 


d(1)"^l0».CH,<1),  +  tC1)1 

while  we  assumed 


if  V3  <  i  <  n+2 

We  observe  that  F(w)  passes  through  points   (0,  0),    (x] ,  p^,    (x-,  p2) 

and  possesses  a  jump  from  (t^,  pn_-,)   to  (**_•,,   1)    (as  allowed  for  an  IFR 
distribution  in  [6,   22]). 


As  it  was  stated  and  proved  form  model  A,  Eq.  (.2.5.4),  objective 
function  L  (2.4.2),  given  feasible  points  (t]. ,  p^),  i  =  1,  2,  ...,  n-1,  is 
indeed  maximized  by  function  F(w)  defined  above  by  Eqs.  (2.5.39-40)  which  1 
a  set  of  the  transformed  exponential  distributions. 

The  feasible  region  of  F(t)  v.s.t.  diagram  for  the  set  of  IFR  distri- 
butions which  pass  through  points  (0,0),  (x.| ,  p-|),  and  (x2>  p2)  is  shown  in 
Fig.  6.  Between  t  =  0  to  t  =  x,  the  region  is  the  area  confined  from  above 
by  F,(t)  •  1  -  e  ]  ,  for  the  same  reason  stated  for  Model  A.  The  region 
is  confined  from  below  by  Max{F2(t),  0.},  where  F2(t)  is  the  t.e.d.  con- 
necting ponts  (x.|,  p^  and  (x2,  p2)  and  according  to  Eq.  (2.5.5)  is  given 
by 


_,Ct-Xl) 

-  O-Pi)  e 


-a2^"V 


The  reason  for  this  is  that  any  point  (t.,  pi )  below  F2(t)  will   be 
connected  to  point   (x, ,   p, )   by  an  IFR  distribution  which  crosses  t.e.d. 
F2(t)  from  below.     But  according  to  property  1    (2.3)  no  IFR  distribution 
can  cross  a  t.e.d.   at  two  points    (e.g.    (x-, ,  p-,),    (x2,   p2))  first  from  below. 
From  left  and  right  the  region  is  confined  by  t  =  0  and  t  =  x,   respectively. 
Between  t  >  x^and  t  *  x2>   the  region  is  confined  from  above  by  t.e.d.    F2(t) 
since  according  to  property  2-a   (2.3),  an  IFR  passes  through  three  points 
(e.g.    (xr   Pl),    (t.,   p,,),    (x2,   p2))  only  if  (t. ,   p.)   lies  on  or  below  the 
t.e.d.   through   (x] ,   p] )  and   (x2,  p2),   i.e.,   F2(t).     The  region  is  confined 
from  below  by  t.e.d.   F^t)  which  connects  points   (0,  0)   to  (x1 ,  p^.     The 


-a,t 

reason  that  {t-,   Pi )  cannot  be  below  F^t)  =  1  -  e    is  that  any  IFR  dis- 
tribution which  crosses  at  point  (t  =  x1 ,  p^  should  have  a  higher  failure 
rate  at  this  point  than  1  -  e    .  So  no  point  Ctj>  Pi )  for  t  >  X|  can 
lie  below  F, (t).  The  region  is  confined  from  left  and  right  by  t  =  x1  and 
t  ■  x2  respectively.  Between  t  >  x2  and  t  =  T  the  region  is  confined  from 
above  by  F(t)  =  1  as  it  is  obviously  clear  and  from  below  by  Fg(t) .  The 
reason  that  Ct,,  p.j )  cannot  be  below  F2(t)  is  that  any  IFR  distribution 

which  crosses  at  point  (x,,,  p2)  should  have  a  higher  failure  rate  at  this 

-a2A 
point  than  1  -  e    .  So  any  point  [ti,   p.)  for  t^  >  x2  cannot  lie  below 

t.e.d.  F«(t).  The  region  is  confined  from  left  and  right  by  t  =  x.,  and 

t  =  T  respectively.  In  addition  to  the  above  limitations  the  point  (t^,  p..), 

must  satisfy  the  equations  given  in  (2.5.39) 

Maximization  for  given  policy  can  now  be  carried  out  by  dynamic  pro- 
gramming methodology.  The  search  extends  over  all  p.  =  F(t.),  i  =  1,  2,  ...,n 
As  for  Model  A,  we  define  Ji ,  as  loss  or  maintenance  cost  on  the  interval 
between  (i-l)th  and  (i)th  inspection  [t-n,  t.],  which  is  also  a  function 
of  the  failure  probability  at  (i-l)th  and  ith  inspection.  That  is 

Ji  -  Ji  (t._},   Pi_r  tr  p.)  =  [I  +  a(ti  -  t.^)]  F*(ti._-| ) 


If  J*  =  Max  J.,  we  obtain  from  Eqs.  (2.5.39)  that 
1   F(w)  1 

J*  -  J*(t,  ,,  p  ,,  t,,  p.)  -  Max  J.(t.  ,,  p  ,,  t 


Feasible  Region 


Fig.  6.  The  Feasible  Region  for  Points  (t-.p^  for 
Model  B. 


F*Ctf]  -  F*(t.    ,] 
|D  +  at^   -  t^,)]  f*(tiA)  +  a 


if  i  *  1,2,   ...,  ik,  m-,+2,   ....  iBg,  m2+2,...,n- 
(1-p,)   -   F*(t     )       F*(tm  +1)   -    (1-p,), 


■  [I  +  a(t1  -  t..1)]F*(ti.1)  ♦  a(- 


if  i  »  nu+1, 
,CH»,)-F*lO       __, 


=  [I  ♦  a(t     +1   -  tra  )]  F*(tm  )  * 


F*(V)"°'P2)i 

i if  i   =  m,+l   and  m, 


0-P2)-F*Ct     )       F*(t         )   -   (l-p2) 
■  [I  +  a(tn.   -  t^,)]^.,)  +  a  ^— *-+  \ 1 


it   i   -  mg-r,    dim  m2  t  my 

=  [I  +  a(ti   -  t^)]  F*^.,)  if  i   ■  n. 

As  mentioned  for  Model   A,   the  first  expressions  on  the  right  side  of 

equations  in   (2.5.44),   i.e.,   [I  +  a(ti   -  t^)]  F*(t._1)  are  obviously  the 

share  of  J.   from  the  first  expression  in  the  expected  total  maintenance  cost 

L  given  by  Eq.    (2.4.2),   i.e.,   I  "[     F*(t.)  +  a"j     F*(tJ(t     ,   -  t.).     The 
1-0  1  i=0 

second  expression  is,  a  /   F*(w)  dw,  which  is  also  the  share  of  J,  from  L, 


t  -T 

a  /   F*(w)  dw.  For  values  of  i,  except  for  i  =  m,+l  when  m,=0,  the  cal- 

0  t.  l  Z 

culation  of  the  values  of  a  /   F*(w)  dw  is  the  same  as  for  model  A.  But 

for  i  =  m,+l  when  t.  >  x„  we  have 

S  xi  \  s 

I       F*(w)  dw  =  /   F*(w)  dw  +  /  F*(w)  dw  +  /  F*(w)  dw.       (2.5.45) 
*1-1         Vl         xl  x2 

Substituting  the  values  of  F*(w)   from  Eqs.    (2.5.38-40)   into  equation 
(2.5.45)  we  have 

dw  + 


/     (1-p,)  e 


i    ~c(m  +31    lw   "   a(m  +1V 
J    e     (  2     J  (V3) 


■C(m1+l)Cw  "  dCra1+l))     ] 


(1-Pt) 


■<W3)  vw  "  >,+3)J 


"tra^U^i-l   *  d(m1+l))         -CCm1+l)(xl   "  %,+!)> 


"c(m2+3)(x2  "  d(m2+3))         -c(m2+3)    (ti   "  d(m2+3)> 
C(m2+3) 

^I'^Vl)1       O-PJ-O-PJ 


c(nyi) 


C(m2+3) 

since  transformed  exponential   distribution  F(w),  which  connects   (p-,,  x, ) 
and  (p2,  x2)  has  already  been  given  by  Fgtt)  as 

F2(t)  -  1   -   (1   -  P])  e'*2     "Xl     , 

and  also  according  to  recodification  (2.5.38)  and  (2.5.40)  we  have 


F*Ctfm   ,)   =   F*(t) 

F*(t,m  +n)  •  F*(x   )  =   (1-p .) 

c(Vl)  =  c* 
F*(tm   )   -    (1-p,)       „  .„  (1-p.)   -   F*(t      ..) 


In  Model   B  also  we  define  K.    (t.,   p.,   c.  +  ,)  as  minimax  expected 
:enance  cost  during  time  interval    [0,   t,]  with  F(tj)  «  pi  when  the 
leter  of  the  t.e.d.   connecting  (t.,   p.)   to   (t-+1>   P-+1)   is  less  than 

R,    (tr   pr   c2)   -J*(t0.   p0.   tr  p,)  =  J*  CO,  0,   tr  P]) 
K^V  Pr  c.+1)  >  K1,l    (.t._r   p._y  Cl)  +  J*  (t._r   p._r  t., 


■jS.fl'v, 


Let  Kn(T)  signify  the  loss  for  the  entire  time  period  T  at  F(T)  =  1  and 
t.e.d.  with  parameter  c  .,  ■  ■ 

Kn^  =  Kn  ^n^'Pn^'VT  ") 

'   Kn-lCtn-T  Pn-T  cn>  +  Jn{tn-T  Pn-1'  tn'  pn} 

•  L.  (2.5.47) 

In  this  model  as  for  Model  A,  we  have  three  state  variable,  i.e.,  t., 

p^  and  z.+y     There  are  two  control  variable,  i.e.,  t.  ,,  p._,.  Minimaxation 

will  be  carried  out  in  three  phases.  In  the  first  phase  the  minimax  expected 

maintenance  cost  fromtimetg  =  0  to  the  time  of  first  inspection  in  the 

interval  [x, ,  x?),  i.e.,  K   ,  (t   .,  p   ,,  c   ,)  and  the  minimax  expected 

i  ml  +  l   V   V   mlV 

maintenance  cost  from  time  tQ  =  0  to  the  time  of  first  inspection  in  the 

interval  [x„,  T],  i.e.,  K  .,  (t   , ,  p   .,  c  .,)  will  be  obtained. 
c  T1   V   V   m2  Z 

In  the  second  phase  K  ,,  (t   , ,  o     .,,  c  ,,),  i.e.,  the  minimax  loss 
m2  i   m2-M   nrigt-i   m2+^ 

from  tQ  =  0  to  the  time  of  the  first  inspection  t  t]  >x.  will  be  optimzied 

for  all  possible  values  of  t   ,,  p   ,  and  c  , ,  =  c  ,,.  In  the  third  and 
m2+l   m2+i     m2+^    m2+l 

final-  phase  Kn  (tn=m  +.  =  T,  pn=m  +.  «  1,  cp+1  =  «)  =  L  will  be  minimaxed. 
In  the  second  phase- the  minimaxation  will  be  achieved,  taking  the  minimax 

K   ,  values  for  state  variables  at  time  t  =  t  ..  as  the  K  values  for  the 
ml  '  "l  ' 

first  stage  of  the  second  phase  and  in  third  phase  also  in  the  same  way  the 

minimax  K  .,  values  for  state  variables  at  time  t  =  t  , ,  are  taken  as  the 
m2+l  m2+l 

K  values  for  the  first  stage  of  the  third  phase.  The  value  of  L  will  be 
calculated  for  all  possible  number  of  stages  in  phase  1,  phase  two 
and  phase  3  and  the  optimum  value  of  L  will  be  found  by  a  simple  search 
among  these  values. 


First  optimization  phase  -  Let  the  minimax  optimized  loss  for  time 
period  (0,  t.~\  at  F(.t.)=pi  when  the  parameter  of  the  t.e.d.  connecting 
(t.,  p..)  to  (t.+r  p.  +  1)  is  equal  to  a  certain  amount  c.  +  „be 

Ki(V  pi'  S'  +  l'  =  Min  Max  Ki  (ti'  V  ci'  (2.5.48) 

where  P  is  the  inspection  policy  (2.2.1)  and  F  is  failure  distribution  given 
by  Eqs.  (2.5.39).  Now  by  using  dynamic  programming  procedure  we  obtain 

K*  (t0-0,  P0=0,  cQ+1)  =  0,  (2.5.49) 

K?  (tr  Pr  c2)  -  J*  (0,0,t1)Pl)  .  I  ♦  at,  ♦  a  F*(tl]  "  1  ,  c*  <  ». 


P"(t-, )  <  p1  <  p'Ct,) 

K*(t2,  p2>  c3)  -  Min  Max{K|(tr  P] ,  c2)  +  J|(tr  P] ,  tj,  P2)} 

0<t2<Xl  (k^-rtj       v(t2,   p2,   c3,   t,)  <  p,   <  G(t2,   p2,   t]) 

p"(t2)  <  p,  £p'(t2) 

Ki  (V  pi'  <W  -  Min  M«tK?.i(t1.l.P1.l.e1)  +  J*(t._1,Pi_1,t..P.). 

0<ti<x1  0<t1.1<t1     v(tl  >P1  ,c1+l  .t^,  )<Pi_1<.G(ti.,pi,ti_1 ) 

P"(ti)<pi<p'(ti) 
K^+1    (tmi+1,Pmi+1,cmi+1) 


in  Max{K*  (t     ,pm  ,c*)  +  J*  +1(t     ,pm  ,t     ., ,pm  ..)} 

H     ml     ml     '  ml    '    ml     ml     V     V 

%^         v^^.p^^,^)^^^) 

n>,  =0,1,2,    ...,  N., 


LafcM 


Km  +1  Ct_  +1  ,    Pm  +1 ,    c  ) 

m2   '     V1       V       m2    ' 

x2-tm2+l  1  T 


Max{K*  (t     ,p     ,c*)  +  J*     .(t     ,pm  ,t     ..,   pm  ..)} 
in-,      in-,      m,      i  hut  i     m       ra       iiiTi       iriiTi 


'2    '     ml     ml     ml    '       V 


'X   *1 


'"1 

m1  -  0,1,2,  ...,  N., 
whsre  N^  is  the  maximum  possible  number  of  inspections  in  the  interval 
CO.x,)  andCl  <c2<  ...  <c*<cvl,   ...  <<^<cg<c^+,<  ...  <cn., 

since  the  F(w)   distribution  is   IFR.     Also  because  at  point   (Xj.p,),   F(w) 
crosses  F(u)  =  1   -  e     1     from  below,  we  should  have  o-j   <  c?  <  c     +, 
P'(t.)  and  P'^t,)  are  defined  the  same  way  as  in  Model   A  (See  Fig.   5.). 
But  their  values  are  given  by  (See  Fig.   6.) 


p'(t,) 


MaxiO,   1   -   (1-Pl)  i 


■A-V, 


id  at  ti   =  T  we  have  p"(T)  =  1. 

The  definition  of  G  and  v  functions   for  Eqs .    (2.5.49)  are  the  s 
;  for  Model   A  and  with  a  similar  reasoning  we  have  (See  Fig.   5-6): 


S(t1,    Pv    tM)    =    1.    -    (1    -    P^    ! 


i.PiXt,  -  t.  .) 


Loggd-p,-) 


GQ(tml)   is   the  particular  value  of  G  fur 
je  is   (See  Fig.   5-6). 


"i^-rlt-T  or0}; 

0  <  t,   <  x,  (2.5.55) 


v^m  +i'Pm  +Ttm  )  is  the  particular  value  of  v  function  at  t 


v(tm  +1,pm  +1,t     )  =  Max{l    -   (1   -  p,)e       '  '  ,0}.         (2.5.57) 

The  value  of  Ci,   i.e.,  the  parameter  of  the  t.e.d.   connecting  point 
(t|,   p^   to   [tj_y   P^).  according  to  Eq.    (2.5.5)    is  given  by 


i  the  value  of  c    +,   is  given  by 


"v        y+i  ■ 


We  now  choose  for  each   (t    +, ,  pm  +,)  that  in,   =  mf  which  renders 


cost  from  til 

Kyi*  ■"»**" 

before  time  t 
Second  optimi; 

expected  main 

second  phase  of  the  optimization. 

K*J+?(tm*+2,  p_*+,,  c  ,+,) 
ml   2    ml  Z      ml  -      V3 

xl   <  W+2  <  x2 


•  P  *+t   c  *     )  is  the  minimax  expected  total  maintenan 

<Y'       Vl 
.  =  0  to  t  *+1  where  tm*+1   is  larger  or  equal   to  x1  wit 
.     Also  m?  is  the  optimum  number  of  checks  or  inspectio 


on  phase  -  The  iterative  procedure  continues  taking  the 
each  .point  (tm*+1>  Pm*+1)»x1  <  tm»+1  <  Xj,  as  the  miniir 
ince  costs  up  to  that  point  for  the  first  stage  of  the 


=  Min  Max{K-+1<y+1,   p^,  c^J 

^mf+l'W     ^+2'   pm*+2'   cm*+3'   ^)<"«fH<8{ti^«''^*2't*fH) 

u™>2(Vr  Vl'  V2'  V+2)} 


K„;+i(V+i,  p_f+i,  c       ) 
mrJ  mi J   nii J    rarJ+1 

P"(tmrJ)<Pmfj<P'(tmr+j) 

=  Min  MaitfISS**j-l<*mT+J-T   f 

Vt  <t  v(t  p  c 


«<W  Pmf+j.  u,l 


(2.5.62) 
i        m2+l 

=  Min  fo^K^m   'Pm    'C?5    +  K  +1(t_    .P.    .t     +1,    pm   .,)} 

m2     m2     mz     c  m2    i      m2     m2     m2    I        m2    i 

I         m2         2  m2+l       m2+l       m2+l       m2  m2         0     m2 

"2"1*"?'   2+mt'    ••"N2  +  "1?' 


N2  is   the  maximum  possible  number  of  inspections   in  the  interval 
[xv  x2]  and  cm<r+1   <  cM»+2,    ....  <  C    +1   since  failure  distribution  F(w) 
is  IFR.     p'ftj)  andp"^)   are  given  by  (2.5.2.13)  and  (2.5.2.14)  re- 
spectively.    The  value  of  G  function  is  given  by  (See  Fig.  5.6) 

-c(t       .,   Pm*+i)(V+1.l   -  x) 

(2.5.63) 

c^j-  vj'  =  *  v+j  -  xi —    '  (2's,64) 


GQ(tm  )   is   the  particular  value  of  G  function  at  point  t       (See  Fig. 
5.6.) 

e  Z  .  (2.5.65) 

The  value  of  v  function  is  given  by 

y+j+i^+rWj-i) 


Al        Lm*+j       V 
v(tm  +1,   pm  +1,  t     )   is  the  particular  value  of  v  function  at  t 


cm  +l(x2"tm  >  "Vm 

v(tm  +T   pm  +1»tm  >  =  Max{1   "   (1    "  P2>  e     2  2       or  1   -  e  2 

V       m2+1     m2  2 


where  1    -  p 

^eC-T^pf) 


The  value  of  cm»+.,   i.e.,   the  parameter  of  t.e.d.  connecting  point 
(tm*+j'  P«ftjJ   t0   (tm*+j-l>  y+j-l)'  ^cording  to  Eq.    (2.5.5)   is  given 


Hfi-J-I 


Now  we  choose  for  each  (tm  +1 ,  p      ^  that  m2  =  m*  =  m*  +  j*  wi 

renders   K**       a  minimum.     That  is 
V 

^(y,.  ty.  ^)  -  H1n{^+1(t^+1.  PV1,  Vl) 
X2-V1<T  ^  -  1  +-f.  a  +«f m2  -• 

m2  -  off  c-  '  - 

^+l.n(T,l,-)-M1n{K-+1(T,l,-) 

m2=l+m!,2+m!,...,N2+mf 

or       ^(T,   1,.)}    , 

m2  -  mi   =  0,1,2,...,   N1  (2.5.71) 

where  Km  +1    is  calculated  in  the  first  optimization  phase  for  m2  -  m-,   or 
equivalently  when  no  inspection   is  performed  in  the  interval   [x, ,  x,). 
Kmf+1    (tm*+T   pm*+l'   cm*+l }   1s   the  minimax  expected  total  maintenance  cost 
in  the  interval    (0,tm,+1]  where  T  >  tj])*+1   >  x2  with  F(tm,+1)  -  pm*+]. 
Also  m?j  ■  m^+j*  where  m^  is   the  optimal    total  number  of  inspections  at 
the  end  of  second  stage  and  j*  is  the  optimal   number  of  inspections  only 
in  the  second  stage. 


Third  optimization  phase  -  The  iterative  procedure  continues  taking  the 
Km*Ii  values  at  each  point  (t  tll,  p^.J.x,  <  t*,,   <  T  as  the  ir 
expected  total  maintenar 
the  third  phase  of  optin 


K*J*2  (t  *+2,  p       _,  c       ,) 
m2  d       m2l       m2  2       m2  3 

x,  <  t  ^.  <  T 


mZl  mZZ  m2+2 

•  Min  MaxfK^jMt  *+1,Pm*+1,   cm*+1 ) 

m2    '     V     m2   '       V 

X2-V<1:iti*+2       v(tm*+2'    pm*+2«    Cm*+3'    tm*+1  )<Pm*+l<G(tm*+?'    Pm*+?>    V+l) 
<:     mz     mfd  mfd       mfd       m2+j       nyi       m2+l         m2+2       m2+2      m2+l 

+  Jm*+2   'V+T   Pm*+T    W'   pm*+?^}' 
m2^      nyi       m2+l       m2+^       m^l 


m*+z       m*+i"     m*+2.'     m|+2.+l 

x2  <  W  *  T 

1  m2+«. 

P"(t_*+J   <  pm,.,   <  P'(tm*,j 
m2  m2   1  ryi 

=  Min  Max{KB?*i-l(tm*n-T  pm*+£-r  W' 

™2;   '     m2  z   '       m2  2   '       V* 

2  m|H-l<  m*n  m*+i'   m*+i'   m*+i+l'   m*n-l    <pm*+i-l 


+  J**+„(t ..,  ..p..,  ,,t m*+0,pm*+ .)} 

2  2  L  '     m2  *  '     m2  *     m2  z 

£  =3,4,...,  n-m*-l.  (2.5.72) 

In  particular,   the  term  for  the  total  optimal    loss  with  n  -  m*+s. 

<£**+.   (T.I.-) 

n  m2  * 

=  Min  Max{ 

x^W.  i<T     *(***».,  l'  <  Pm*+«  i  <  G  (*„,*+,  i) 
<:     m2+4-l  m2  t  m2  i   i  mpi-i 

+  J*  (tm*+9   ,,   pm*+.   ,,t  ...    =T,pn_m)k+  =1)} 

n=m*+i      m2  *■   '       ra2  l'[     n'm2  l        n'mZ  l 

i  =  1,3,... ,N3,  (2.5.73) 

where  N3  is  the  maximum  possible  number  of  inspections  and  G  and  v  functions 

are  given  by  (See  Fig.  5-6.) 

-c(V+,,Pn]lt+£)(tm,+z_rx2) 

sWWVm  i>  =  ]  -  (]-p?)e      2      2       2 

mZ*    mZl    m2 nl  2 

(2.5.74) 


2  2  VJ"X2 


id  G(t  *.     J  is  particular  value  of  G  funct 
;  given  by 


v(tm*+£,pm*n'cm*n+l  'V+n-l' 


'"2  x     m2  *     '"2  2 


■  Max{l   -   (1   -  p  »     )  e     2  2  2  (2.5.76) 


1    -     (1     -    P-)    6 


1  \ 


»(w»  ■  i  -  d-p2) e 


The  value  of  c  *,.  can  be  found  according  to  Eq.   (2.5.5)  for 
cn  t  cm*+l  by 


subject  to  a,  <  c  tll   <  c  ».-  <   ...   <  c     !__*..    -,   since  the  failure  dis- 
2        m2+l  -    m2+2  -  -    n-1-iy-i-l 

tribution  should  be  IFR  and  cm<:+1    is  given  by  Eq.    (2.5.68)   for  m2=m*. 

The  optimal   loss  L**  at  n*  =  mi  +  i*  =  m*  +  j*  +  £*  where  n*  is  the 
optimal   total   number  of  inspection  in  the  interval    [0,T]  and  i*  is  the 
optimum  number  of  inspections  in  the  interval   [x-.T]  is  then  found  from 

►*(T,1,»)  or     K(T,1.-)         ] 


;  given  by  Eq .    (2.5.71). 


The  P**  and  failure  distribution  F**  defined  by  Eqs.    (2.1.3.4)   can  also 
be  found  if  needed  by  a  backward  recursive  procedure  from  the  computer 
print  out. 
2.6     DISCUSSION 

In  this  chapter  the  mathematical  meaning  of  a  minimax  policy  was 
stated  followed  by  the  statement  of  the  present  problem  which  consists  of 
two  models  A  and  B  according  to  having  information  about  one  point  or 
two  points  of  an  increasing  failure  rate  (IFR)   distribution  of  a  system 
respectively.     The  special    properties  of  IFR  distribution  utilized  in  the 
logic  of  the  formulation  of  the  model  A  and  B  were  proved  and  explained 
by  property  No  1   and  property  No  2. a  and  2.b.     The  form  of  the  objective 
function  suitable  for  a  recursive  relationship  was  presented  and  derived. 
Both  model  A  and  model   B  have  been  formulated  by  functional  equations  of 
dynamic  programming.     Three  state  variables,   i.e.,  the  last  inspection 
time  t-,   the  cumulative  failure  probability  of  the  system  F(t|)  ■  p^  up  to 
time  t^   and  the  parameter  of  a  transformed  failure  distribution  passing 
through  points    ( t^ , pi )  and  (t<+,  ,p.+, )  or  c.+1    represent  the  state  of  the 
system.     The  control  variables  are  t.   ,,  the  timing  of  the  previous  in- 
spection and  p^-j,  the  failure  probability  of  the  system  at  t._-|.     The 
stages  represent  the  number  of  inspections  which  ranges   from  2  to  maximum 
possible  number  of  inspections   in  the  interval    [0,T].     The  minimax  policy 
first  maximizes  at  each  stage  and  for  a  fixed  state  vector  (t- ,p^,c+, ), 
the  total  expected  loss  from  initial   state  vector  (t0=0,pn=0,cn+, )  to 
the  final   state   (t^.P^.c^,)    by  the  choice  of  the  control   variable  p.    , 
and  then  minimizes  the  total  expected  loss  by  the  choice  of  control 
variable  t.    j .     As  it  is  clear  the  number  of  stages  should  be  also  optimized. 


This  has  been  done  in  model  A  in  two  phases,  i.e.,  in  phase  one  the  optimal 
number  of  inspections  or  stages  in  the  interval,  [0,x,)  is  obtained  and 
in  the  second  and  final  phase  the  total  optimum  number  of  inspections  and 
their  timing  is  obtained.  In  model  B  optimization  has  been  accomplished 
in  three  phases  where  the  optimal  number  of  inspections  in  the  intervals 
[O.X,),  [xrx2),  [x2,T]  are  obtained. 

In  Chapter  3  the  computational  details  and  procedures  will  be 
presented. 


CHAPTER  3 
COMPUTATIONAL  PROCEDURE 

3.1  INTRODUCTION 

In  this  chapter  the  numerical   procedure  for  the  computer  program  is 
presented  for  each  of  the  two  models,   i.e.,  Model   A  and  Model   B.     For 
each  model   first  the  range  of  the  possible  values  of  state  and  control 
variables   is  stated  and  then  the  computational  and  numerical   procedures  is 
followed  step  by  step  with  reference  to  the  computer  program  in  Appendix  B. 

3.2  MODEL  A 

In  Chapter  2,   there  are  three  state  variables  and  two  control   variable: 
in  the  formulation  of  the  problem.     The  state  variables  are: 

1.  The  inspection  or  checking  time  of  the  ith  inspection  or  checking 
denoted  by  t. . 

2.  The  cumulative  failure  probability  of  the  system  at  ith  inspection 
denoted  by  p. . 

3.  The  parameter  of  the  transformed  exponential   failure  distribution 
connecting  point  (tj,  pi )   to  point  (ti+],  pi+1)  on  the  F  v.s.   t 
diagram  denoted  by  c+,. 

1.  The  inspection  time  of  the   (i-l)th  inspection  denoted  by  t.    ,. 

2.  The  cumulative  failure  probability  of  the  system  at  (i-l)th 
inspection  denoted  by  p.,. 

ti   and  tj_i    according  to  the  assumption  f  in  (2.2)   are  discrete  values 
ranging  from  0  to  T  the  maximum  life  time  of  the  system,     p.   and  p. 
haye  continuous  values  whose  range  is  given  by  p'  and  p"  functions  defined 


i  the  range  of  its  values  is  giv 


LCi+1<R2     , 

;  of  R,   and  R.  are  given  by 


1-Pi 
l°9e<HT> 
t.-x 

L°3e(& 


Loge(y; 


in"2  = 


(3.2.3) 
the  properties 


The  reason  for  the  above  values  for  R1  and  R~  is  based 
of  exponential  and  IFR  failure  distributions  explained  in  (2.3)  and  on  tl 
same  line  of  reasoning  used  in  (2.5.2)  for  stablishing  the  feasible  regit 
for  IFR  distributions  passing  through  points  (0,0),  (x^p^  and  (x-.p,)  < 
the  F  v.s.   t  diagram  (See  Fig.  6)  with  the  difference  that  here  the  feas 


region  is  given  for  the  IFR  disl 
(t^.Pj)  and  (x,p)  for  tj  in  the 
(t^.pj   and  (T,l)   for  t.    in  the 


tions  which  pass  through  points   (0,0), 
interval    (0,x)   and  through  points   (x,p) 


The  input  vector  elements  are: 

a  ■  The  cost  of  undetected  failure  per  unit  time. 

I  *  The  cost  of  every  inspection  with  the  exception  of  the  inspection 

cost  at  t  =  0  where  1(0)  =  0. 
x  =  the  time  at  which  the  cumulative  failure  probability  of  the 

system  is  known, 
p  =  The  known  cumulative  failure  probability  of  the  system  at  time  x 

r  «  {t0  «  o,  tr  t2 x,  ...  tN  =  n 

"  The  set  of  possible  inspection  times. 
Pit,*  Number  of  increments  of  pi    &  pi_]   in  the  feasible 

region  at  t-    &  ti_]    respectively 
AL^  Number  of  increments  of  ci+1    in  the  feasible  range 
given  by  Eqs.    (3.1.1-3). 
The  computational  procedure  consist  of  the  following  steps; 
Phase  1: 

1.  Divide  the  feasible  range  of  p.   values,   given  by  p'   and  p"   functions 
in  Eqs.    (2.5.15-18),  at  each   time  t,  where  t-    is  an  element  of  the  set 
T',   into  PN,    increments. 

2.  Divide  the  feasible  range  of  cf+1    values,   given  by  R]   and  R2  functions 
in  Eqs.    (3.2.1-3),  at  each  point  (t^.p^ )  on  the  F  v.s.   t  diagram,   in  to 

3.  Assume  no  inspection  in  the  time  interval    (0,T)  with  only  one  inspection 
at  time  T,   i.e.,  set  the  stage  number  m  equal    to  0. 

4.  Calculate  the  maximum  expected  total  maintenance  cost  when  m  inspections 
are  perfromed  in  the  interval    (0,x),   i.e.,   K*=  +,(T,1,=)   and 


Km+l(Vl'pm+l,cm+l)  for  a11  Possible  values  of  Wpm+1  and  cm+l  from 
Eq.    (2.5.7)   for  i   =  m+1   and  set  of  aquations   (2.5.14). 

5.  Set  current  m  =  old  m+1   and  calculate  Kt(t. ,p, ,c.+, )  given  by  Eq. 
(2.5.14)  when  i    =  m  and  for  all   possible  values  of  tj ,  pi   and  ci+,. 

6.  Repeat  step  4  and  compare  the  R*^|,»)  and  K*+]  (Vl  .p^.c^  )   values 
for  the  old  m  values  with  their  values  for  current  m  respectively  and  choose 
that  value  of  K*=m+1  and  K^foreach  point  (tm>p J  which  is  smaller  and 
discard  the  larger  values. 

7.  Compare  the  maximum  value  of  m,   i.e.,   N,  with  the  current  m.     If 
N-|   >  m+1,  go  back  to  step  5  and  if  N-,    <  m+1   go  to  step  8. 

8.  Take  that  value  of  R*^ ,-)  and  K^,   pm+] ,   z^)  at  each  point 
^m+l'Pm+l)  which  1's  ml'nimum  for  all   values  of  m  =  0,   1,    ...,  N,   and  denote 
them  as  i<n=m«+1  (T,l ,»)  and  K**+]  ( tm«+1  ,  pm<:+1 ,  C^,)   respectively,     m* 

is  the  optimum  number  of  inspections   in  the  time  interval    (0,   x). 
Phase  2 

9.  Consider  the  K*J1(t^+1>   pm<+] ,   C|n,+1)   values  as   the  minimax  expected 
maintenance  cost  upto  the  time  of  the  first  inspection  in  the  time  interval 
[x,T)  or  the  first  stage  j  =  1   of  the  second  phase. 

10.  Calculate  K**m*+j (T,l ,»)   from  Eq.    (2.5.32)   for  j  =  1 . 

11.  Set  current  j  *  Old  j+1 .  Compare  the  maximum  possible  number  of  in- 
spections in  the  interval  [x,T),  i.e.,  N2  with  the  current  j.  If  N-  >  j 
go  to  step  12  and  if  N2  <  j  go  to  step  14. 

12.  Calculate  K*J+j(tm,+j,  p^+j,   c^+j+1)   from  Eq.    (2.5.31)   for  all 
feasible  values  of  t  *+ . ,   p  ,+  .  and  c  ,    .   ,. 

13.  Calculate  K**m*+J-(T,1 ,-)  from  Eq.  (2.5.32).  Compare  K**  t+.  values 
for  old  j  values  with  its  value  for  current  j  and  keep  the  smaller  value 
and  discard  the  larger  values.     Then  go  back  to  step  11. 


14.     Find  L**,the  upper  bound  for  the  optimum  expected  total  maintenance 
cost  with  known  system  failure  distribution  from  Eq.    (2.5.37). 

At  each  stage   i  and  each  state  (t^ ,  p^,   C,+-.)   the  values  of  the 
ti-l'  Pi-i'   1-e-'   the  control   variables  of  stage  i  which  minimaxes 
Mtf,   Pi ,   ci+,)  are  recorded  and  printed  out.     The  value  of  Cj    is   taken 
as  the  first  larger  value  of  R,  +  a.r  where  r  is  the  number  of  increments 
and  a  is  the  size  of  the  increment  obtained  from 

R2-R 
A'-fq1.  (3-2-4) 

respect  to  cc  given  according  to  Eq.    (2.5.5)  by 
1-P,- 
CC  -  t     .  t     1'1      ■  (3.2.5) 

The  m*  values  for  the  first  stage  of  the  second  phase  are  being  recorded 
and  printed  out  at  the  beginning  of  the  second  phase,  i.e.,  when  j  =  1. 
Finally  the  optimum  total   number  of  stages  or  n*  ■  m*+j*  is  recorded 
together  with  L**.     The  policy  p**  and  failure  distribution  F**  defined 
by  Eqs.    (2.1.3-4)   can  then  be  found  easily  if  needed  by  searching  backward 
from  stage  n*  =  m*  +  j*  to  n*-l,   n*-2,    ...,   until   m*+l .     Then  the  search 
continues  from  stage  m*  to  m*-l,  m*-2,   ...,  until   the  first  stage. 

3.3     MODEL   B 

In  this  model   also  there  are  three  state  variables,   i.e.,   t- ,   p..,  c.    , 


and  two  control  variables,  i.e.,  t.  ^,  p^_^  which  are  defined  in  the  same 
way  as  for  Model  A.  Here  also  t|  and  t,  ■.  are  discrete  values  according 
to  assumption  f  in   (2.2)  and  range  from  0  to  T  the  maximum  life     time  of 


the  system,  p^  and  p.,  have  continuous  values  whose  range  is  given  by 
p'  and  p"  functions  defined  in  Chapter  2.  ci+1  is  also  continuous  and 
its  range  is  given  by 


i  the  values  of  R,   and  R„  are  given  by 


IW  P^   = 


-^ 

■1-D,0 

Log6 

1"P« 

'1^ 

LogjT—4 
e   I   p2 

1      *1 

-  x2 
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The  reason  for  the  above  values  for  R,  and  R2  is  based  on  the  properties 
of  exponential    and  IFR  failure  distributions  and  on  the  same  line  of  reasoning 
utilized  in   (2.5.2)  for  establishing  the  feasible  region  for  IFR  distributions 
passing  through  points   (0,  0)   (x] ,  p-,)  and  (x2>  p2)  on  the  F  v.s.   t  diagram 
(See  Fig.   6.)  with  the  difference  that  here  the  feasible  region  is  given  for 
IFR  distributions  which  pass   through  points   (0,  0),    (ti ,  pi )  and  (xr   p,) 
for  ti   in  the  interval    (0,   x1 ) ,   through  points    (*j,  p-,),   (tj ,  Pi )  and 
(x2,p2)  for  t.   in  the  interval    (x^x*)  and  through  points   (x2,P2)    ( t .  ,p . ) 
and  (T,l)   for  ti   in  the  interval    (x2>T).     The  input  vector  elements  con- 
sists of  a,    I,  PN,   and  AL,   as  defined  in   (3.2)  plus 


=   The   t 


at  whic 


-.hi 


:   fail 


ributior 


p1   and  p2  • 


2  failure  probabilil 


•>*2 VT} 


■  The  set  of  possible  inspection  times. 
The  computational   procedure  consists  of  the  following  steps; 
Phase  1: 

1.  Divide  the  feasible  range  of  p^   values  given  by  p'   and  p"  functions 

in  Eqs.    (2.5.50-51),  at  each  time  t<  where  t,   is  an  element  of  the  set  T' , 
into  PN,    increments. 

2.  Divide  the  feasible  range  of  c,+1   values,  given  by  R-\   and  R2  functions 
in  Eqs.    (3.3.1-3),   at  each  point  (tj.  p,-)   on  the  F  v.s.   t  diagram,  into 
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3.  Assume  no  inspection  in  the   time  interval    (0,T)  with  only  one  in- 
spection at  time  T,   i.e.,   set  the  stage  numbers  m,  *  nu  ■  0. 

4.  Calculate  the  minima*  cost  K*     At     .,,   pm  .,  ,   cm  ..)   for  all   possible 

ml    '     V1       V1       V1 
values  of  t       .,   p  and  cm   .,   from  Eqs .    (2.5.44)   for  1   «  m-,+1,  and  the 

"r     V         mi+1  i 

set  of  equations    (2.5.49).     Also  calculate  the  minimax  cost 

^m  +i  ^m  +i '   pm  +t   cm  +i)    for  a11   possible  values  of  t     ,.,  pm  ..   and 
V     m2  ]       m2    '       m2    '  m2+1       m2  ] 

cm  +1   from  Eq.    (2.5.44)   for  i    =  iy-1 ,  m2  «  m-,   and  the  set  of  equations   (2.5.49). 

5.  Set  current  i^   =  Old  n^+1    and  calculate  K*^,  p.,   c+,)   given  by  Eq. 
(2.5.49)  when  i   *  m-,  and  for  all   possible  values  of  t.  ,  p.   and  c+,. 

6.  Repeat  step  4  and  compare  the  K*  +,(t     +, ,   p     +, ,   c     +, )  and 
or  old  m,   values  with  their  values   for 


point  (t           p 

)  and  (tm  +Tpm  +1]  * 

hich  is  smaller  and  disca 

rd  the  lai 

7.     Compare  the  ma 

ximum  value  of  mr   i. 

e.,  N,  with  the  current  v 

alue  of  m. 

If  N,   >m1+l,  go  b 

ack  to  step  5  and  if 

N-,   <  m1   +  1    go  to  Step  3. 

8.     Take  that  vaiu 

e  of  KjJ     ,   at  each  po 

int  ttm  +1 ,Pm  +1)  which  i 

S  minimum 

for  all    values  of 

d  denote  it  as  K**+1(t  t . 

1 '  pm*+l ' 

cm*+1)  where  m*  is 

the  optimum  number  o 

1         ml 
f  inspections   in  the  time 

1 

(0.x,). 

Phase  2 

9.  Consider  the  K**  .    (t  *., ,  Pm*+1 .  c  *  ,)   values  as   the  minimax  expected 

1    '         l1       V       "l    ' 
maintenance  cost  upto  the  time  of  the  first  inspection  in  the  time  interval 

[x1 ,   x2)  or  the  first  stage  j  =  1   of  the  second  phase. 

10.  Calculate  K**  .(t       ,,   p  cm  . . )   values  for  all   feasible  points 

V1   V1     V1     V1 

(tm  +1,  Pm  +1)  on  the  F  v.s.    t  diagram  from  Eq.    (2.5.44)  setting  1   «  m2+l 
and  m2  =  m*+l   and  from  Eq.    (2.5.62). 


Set  current  bu  ■  old  m-+l .     Compare  the  maximum  possible  number  of 
pections  in  the  interval    [x-|,   x2),   i.e.,  N2  +  m|  with  the  current 
ue  of  m2.      If  N2+m|  >  m2  go  to  step  12.     If  N2  <  m2  go  to  step  15. 

Calculate  K**+j  (tm*+j,  pm*+j,   c^+j )  for  all   feasible  values  of 
+.,  Pm*+;  and  c  *+.  from  Eq.    (2.5.44)  for  i   ■  mf+j  =  nu  and  from  Eq. 
2^5. 61). 

3.  Calculate  C**i(tL  j.i»  Pm  *n  cm  j.i)  values  for  all   feasible  points 

V1     m2   '       V       m2   ' 
t_  +1.  Pm  +1)  from  Eq.    (2.5.44)   setting  i   =  m,+l   and  from  Eq.    (2.5.62). 

V     V  i 

4.  Compare  K**  ,  values  for  old  m0  values  with  K**  ,  values  for  current 

V  2  m2    ' 

m,  and  for  all    feasible  points   (t     .,,   p     ,, )  and  keep  the  smaller  values 
I  m2+l       m2+l 

and  discard  the  larger  ones.     Go  back  to  step  11. 

15.  Take  that  value  of  K**  ,    (t       , ,   p     ,,,  c     .,)  at  each  point  (t     .,, 

T1      V      m2''      m2   '  "r 

pm  ,,)  except  point  (T,l)  which   is  minimum  for  all   values  of  in,  ■  m*+l , 
m2+l  2         1 

m|+2,    ...,  m^+N2  as  found  in  step  14  and  compare  it  with  the  minimum  value 

of  K  +i(tm  *t   Pm  4.1  •  cm  4.1)  as  found  in  step  6.     Denote  the  smaller  value 

m2    '     V       m2   '       m2    ' 
of  them  at  each  point  (t      , ,  pm     ,)  except  point  (T,l)  as  K*J*,   (t  ».., 
m2    '       m2    '  m2   '       m2   ' 

Pm**i  >  cm*-n)  where  m*  is   the  optimum  number  of  inspection  in  the  interval 
m2   '       V  l 

(0,x2). 

Phase  3 

16.  Consider  C;j,(t»,,   Pm*+i >  cm*+i )  values  as  the  minimax  expected 

m2   '     T1       m2   '       *2   ' 
maintenance  cost  upto  the  time  of  the  first  inspection  in  the  interval 

[x,,T)  or  the  first  stage  %  =  1   of  the  third  phase. 

17.  Calculate  K**(T,1,»)   from  Eq.    (2.5.73)   for  t  -  1 

18.  Set  current  l  ■  Old  £  +  1.     Compare  the  maximum  possible  number  of 
inspections  in  the  interval    [x,,   T),   i.e.,   N,  with  the  current  i.     If 
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N3  >  i  go  to  step  19  and  if  N3  <  i  go  to  step  21. 

19.  Calculate  K*$*     (t_*.,,   p  *..,   c  *...•,)  from  Eq.    (2.5.72)  and  Eq.    (2.5.44) 

m2  »      m2  *      ragf*       m2  s.    i 

for  i  =  mi  +  i  and  mi  t*  m?  for  all   feasible  values  of  t_*.i,  Pm*.£  and 
£  £         I  m2+        m2+ 

cmi!+i  and  c  *.*.. 
m2  m2 

20.  Calculate  k**  *+„(T,l ,»)   from  Eq.    (2.5.73)   and  from  Eq.    (2.5.44)   for 

2 
i  =  n.     Compare  k**  *.     values  for  old  l  values  with  its  value  for  current 

m2  * 
l  and  keep  the  smaller  value  and  discard  the  larger  values.     Then  go  back 

to  step  18. 

21.  Find  L**  the  upper  bound  for  the  optimum  expected  total  maintenance 
cost  with  known  system  failure  distribution  passing  through  points   (0,0), 
(x-,,   p,),    (x2,  p2)  and  (T,l)   from  Eq.    (2.5.78)  where  Kn=m*+1(T,1 ,»)   is 
given  by  Eq.    (2.5.71).     The  value  of  I*,   i.e.,   the  optimum  number  of 
inspections  in  the  interval    [xg.   T)   is   then  obtained. 

As  for  Model  A,  at  each  stage  i    and  each  state   (t,  ,  pi ,   c,+,),   the 
values  of  the  t^,  p,   ,,   i.e.,   the  control    variables  of  stage  i  which 
5  Kj(tj.   p^>  ci+1)   are  recorded  and  printed  out.     The  value  of 
ten  as  the  first  larger  value  of  R-j  +  i.r  where  r  is  the  number 
nents  and  a  is  the  size  of  the  increment  obtained  from 


sspect  to  cc  given  according  to  Eq.    (2.5.5)   by 
I-P, 


The  mf  values  for  the  first  stage  of  the  second  phase  is  b 
recorded  and  printed  out  at  the  beginning  of  the  second  phase, 


j  =  l.     The  m^  values   for  the  first  stage  of  the  third  phase  are  also  being 
recorded  and  printed  out  at  the  beginning  of  the  third  phase,  i.e.,  when 
!■  1.     Finally  the  optimum  total   number  of  stages  or  n*  =  mi  +  i*  = 
ircf  +  j*  +  £*  is  recorded  together  with  value  of  c**.     The  policy  p**  and 
failure  distribution  F**  defined  by  Eqs.    (2.1.3-4)   can  then  be  found  easily 
if  needed  by  searching  the  printout  backwardly  from  stage  n*  =  m?+j*  +  l* 
to  n*-l,  n*-2,   ...,  until  m|+j*+l.     Then  the  search  continues  from  stage 
m*  +  j*  to  m*  +  j*-l,  m|  +  j*-2,    ...,  until  mj  +  1 .     Now  the  search  continues 
from  stage  WJ  to  mf  -  1,  nvf  -  2,  . ...  until   the  first  stage. 

In  appendix  B  the  computer  program  together  with  its  discription 
is  presented.     In  Chapter  4  the  result  for  both  Model  A  and  Model   B  are 
presented  and  comparison  is  made  between  the  two  models  on  the  basis  of 
L**  values  obtained  by  each  model. 


CHAPTER  4 
THE  RESULTS  &  APPLICATION 

4.1  INTRODUCTION 

In  this  chapter  the  results  obtained  from  computer  programming  for 
several  input  data  has  been  presented  and  the  application  of  the  results  ha 
been  discussed  through  an  example  problem  utilizing  the  information  obtaine 
from  it.  The  results  are  presented  according  to  the  following  format. 

1.  The  convergence  and  accuracy  of  the  results. 

2.  Evaluation  of  the  value  of  information  about  the  failure  distribution 
of  a  system  through  comparison  of  the  maximum  expected  total  maintenance 
cost  L**,  applying  Model  A  (location  of  one  percentile  known)  and  Model  B 
(location  of  the  two  percentiles  known). 

3.  Comparison  of  the  actual  optimum  expected  total  maintenance  cost  for 
several  IFR  failure  distributions  with  the  maximum  of  the  optimum  total 
expected  maintenance  cost  L**  for  both  Model  A  and  Model  B. 

4.  Evaluation  of  the  value  of  information  about  a  system  failure  distri- 
bution with  regard  to  the  relative  location  of  the  known  points. 

5.  Variation  of  L**  with  the  inspection  Cost  I  for  Model  A  and  Model  B. 

6.  Variation  of  L**  with  the  cost  per  unit  time  of  undetected  failure  a_ 
for  Model  A  and  Model  B. 

7.  Example  problem. 

4.2  CONVERGENCE  AND  ACCURACY 
FOR  MODEL  A 

The  following  input  data 


=  1.4  $/inspectior 


T1  =  {0,0. 1,0. 2, 0.3, 0.4, 0.5, 0.6,0. 7, 0.8, 0.9,1., 1.1, 1.2, 1.3,1. 4}, 
plus  each  of  the  following  sets  of  values: 


were  given  for  the  computer  program  in  Appendix  B.  The  values  of  F**,  P** 
and  maximum  expected  total  maintenance  cost  for  different  numbers  of  in- 
spections in  each  of  the  subsets  of  inspection  times 
Tj  =  {0,.1,.2,.3}  and  T£  =   {.4, .5, .6, .7, .8, .9,1 . ,1 .1 ,1 .2,1 .3} 
are  given  in  Table  1  where  SS-, ,  SS2  and  SS3  refer  to  the  values  of  F**,  P**  a 
maximum  expected  total  cost  v.s.  number  of  inspections  when  the  number  of 
increments  of  state  variables,  c  and  £  are  given  by  sets  S, ,  S.  and  S3 
respectively.  Table  1  also  shows  that  the  minimum  value  with  respect  to 
number  of  inspections  of  the  maximum  with  respect  to  system  failure  dis- 
tributions of  IFR  type  passing  through  point  (x  =  .4,  p  =  0.180),  i.e., 
L**  is  equal  to  10.1  for  M2  =  5.  P**  shows  a  good  convergence  and  accuracy 
since  shifting  from  set  S.,  to  S^  with  higher  number  of  state  increments  has 
not  changed  the  solution  for  P**.  However  for  F**  and  expected  total  cost 
the  table  shows  that  although  shifting  from  S1  to  S2  and  from  S2  to  S3  in- 
creases the  convergence  and  accuracy  but  not  to  more  than  one  digit  of 


accuracy  has  been  obtained, 
state  variables  increases  tl- 
creases  very  rapidly  since  i 


Surely,  inc 
!  accuracy  b 


:   this 


expected  cost,  in  many  cases  the  level  of 
order  obtained  for  this  problem  (i.e,  yjpy 

FOR  MODEL  B 

The  following  input  data 

a  ■  35.  $/Unit  Time, 
I/Inspection, 
Unit  Time, 


sing  the  number  of  increments  of 
the  computational   time  also  in- 
itinuous  state  variable,   i.e., 
i  upper  bound  for  the  total 
:curacy  on  this  upper  bound  of  the 
<  100  «  99%)  can  be  sufficient. 


T'   =   {0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3,1.4}, 
plus  each  of  the  following  sets  of  values: 


-, 


;AL, 


were  given  for  the  computer  program  in  Appendix  B.  The  values  of  F**,  P** 
and  maximum  total  expected  maintenance  cost  for  different  numbers  of  in- 
spections in  each  of  the  subsets  of  inspections  T,'  ■  {0,0.1,0.2,0.3}, 


T£  •  {0.4,0.5,0.6,0.7,0.8}  and  T^  =  {0.9,1.0,1.1,1.2,1.3}  are  given  in 
Table  2  where  SS1 ,  SS2  and  SS3  are  defined  in  the  same  way  as  for  Model  A. 
Table  2  also  shows  that  the  minimum  value  with  respect  to  number  of  in- 
spections of  the  maximum  with  respect  to  system  failure  distributions  of 
IFR  type  passing  through  points  (x]  *  0.4,  p]  =  0.180)  and  (x-  =  0.9, 
P2  »  0.560),  i.e.,  L**  is  equal  to  9.61  for  M3  =  3.  Again  here  P**  shows 


up 


good  convergence  and  ace 
gher  number  of  state  inc 
r*  and  maximum  expected  c 


ice  shifting  from  i rip 
las  not  changed  the  s 
lues  in  Table  2  show 
ained  by  shifting  frc 


set  S, 


3    S,. 


4.3     COMPARISON  OF  MODEL  A  WITH  MODEL  B 

Figure  7  shows   the  variation  of  maximum  expected  cost  with  the  number 
of  inspections.      It  can  easily  be  seen  that  the  inspections  which  are  per- 
formed after  the  time  for  which  cumulative  probability  of  failure  is  known 
are  more  important  and  their  number  determines  the  upper  bound  on  the  expected 
total   optimum  cost  of  the  maintenance.     Table  3  shows   the  value  of  infor- 
mation for  IFR  distributions  which  pass  through  points   (0.4,   0.180)  and 
(0.9,  0.560).     If  both   informations  are  utilized  then  Model   3  gives  the 
values  for  L**  which  is  equal    to  9.61.      If  only  point  (0.4,  0.180)   is  known 
then  Model  A  gives  L**  equal    to  10.10.     Similarly  if  only  point  (0.9,  0.560) 
is  known  then  Model  A  gives  L**  equal   to  9.64.     Table  3  also  shows  that  the 
additional    knowledge  about  point  (0.9,  0.560)  improves   the  upper  bound  cost 

formation  about  one  point  is  available  the  closer  this  point  to  the  final 
point  (T,l)  the  more  improved  upper  bound  value, L**  can  be  obtained.     This 


I  x,  For  IFR  Passing 
[   Through  Pt.  (0.4, 
|    0.180)  and  Pt.  (0.9, 
I    0.560). 
I  ©,  For  IFR  Passing 
I    Through  Pt.  (0.9, 
]    0.560). 

©,  For  IFR  Passing 
i   Through  Pt.  (0.4, 


12         3         4       5 


0       12         3       4         5 
.  Number  of  Inspections  During  1th  Period  3 

Given  by  T. 


Fig.   7.     Variation  of  Upper  Bound  for  the  Expected  Total  Cost 
with  M.  Number  of  Inspections. 


Table  3.     The  Value  of  Information  about  IFR  Distributions  Passing  Through 
Points   (0.4,   0.180)  and   (0.9,  0.560). 

Utilizing  The  Information(s)  L**  L  **-!** 


P2  -  0.560; 
(Using  Model 


p  «  0.560;     x  =  0.9  9.64  4.55 

(Using  Model   A) 
For: 

a  •  35.         $/Unit  Time, 
1-1.4        S/Inspection, 

T'   =  {0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3}    , 
and      AL,   &  PN,       given  by  set  S-. 


is  clear  again  from  Table  3  since  having  the  knowledge  about  point  (0.9,0.560) 
is  4.55%  more  yaluable  than  the  knowledge  about  point  (0.4,   0.180)  which  is 
closer  to  point  (0,0).     In  section  4.5  more  results  are  presented  regarding 
the  location  of  the  known  point  or  points. 

4.4  ACUTAL  OPTIMAL  COST  V.S.  L** 

The  optimal   policy  P.    «  {t-|,   t-,    ..-,   t  }  and  the  optimal    cost  R.  has 
been  found  numerically  by  using  computer  program  of  Appendix  C  for  several 
different  IFR  failure  distributions  Fi  which  consist  of  points   (t.  ,  p.) 
connected  by  transformed  exponential   distributions.     Tables  4  and  5  give 
the  policies  P\   and  optimal    costs  Ri   and  also  the  ratio  of  Ri    to  L**.     It 
can  be  seen  that  in  both  tables  this  ratio  is  less  than  one  as  could  have 
been  expected  since  by  definition  L**  is  an  upper  bound  for  the  optimal   ex- 
pected maintenance  cost  for  known  distributions. 

4.5  VALUE  OF  INFORMATION  V.S.    RELATIVE  LOCATION  OF  THE  KNOWN  POINT(S). 
A  certain  IFR  failure  distribution  FQ  has  been  selected  for  which: 

p(t  »  0.0)   =  0.000       ;       p(t  =  0.2)   =  0.150     ;     p(t  =  0.3)  -  0.216; 

p(t  *  0.4)   =  0.278       ;        p(t  =  0.5)   »  0.535      ;     p(t  =  0.6)   =  0.700. 

The  set  of  possible  inspection  times  T'   is  given  by 

T'   =  {0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8}    . 

Now  for  the  above  failure  distribution  the  values  of  L**  are  plotted 

in  Fig.   8  and  Fig.   9  assuming  that  only  one  point  (x,  px)  or  two  points 

(x,,  p     )   and   (x?,  p     )  of  the  distribution  Fn  is  known.     Figure  8  shows 
i       x]  i       x2  U 

that  the  value  of  L**  improves  as  the  known  location  x  of  the  100. p     th 

percentile  of  the  failure  distribution  moves  toward  the  maximum  life  time 
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Table  4.  (continued) 

V   ■  {0,0.1,0.2,0.3,0.4,0. 

,   t7  ■  0 

.6  Unit 
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.2  Unit 
.4  Unit 

ime; 
line; 
rime; 
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0.8  Uni 
1.0  Uni 
1.3  Uni 
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1.2,1.3} 
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Table  5.  Optimal  Policy  P.  and  Cost  R..  for  Some  IFR  Distributions  F. 
Pass  Through  Point  (0.4,  0.180)  and  Their  Comparison  with  f 
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a  =  35.  $Unit  Time 
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T'=  Defined  in  Table  4. 


=  0.0  Unit 
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-  0.7  Unit 

Time; 
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Time; 

■  1.4  Unit 

Time 

E.  •  0.4  Unit  Time 
E.q  ■  0.9  Unit  Time 
t]4  =  1.3  Unit  Time 


=  20.  S/Unit  Time 
■  0.8  $/Inspectio 


3.2   0.3     0.4    0.5    0.6 
X  =  The  Time 
where  px  =  F(x)  is 
Known  in  Units  of  Time 

-iation  of  L**  with  the  Time  x  a 
=  F(x)  is  Known. 


of  the  system.     The  same  trend  can  be  seen  in  Fig.   9  for  constant  value  of 
x-.     In  addition  to  this,  Fig.  9  shows  that  L**  improves  as  the  location 
of  x2  moves   toward  the  maximum  life  time  of  the  system. 

4.6  VARIATION  OF  L**  WITH  I 

Figure  10  shows  the  variation  of  L**,  the  upper  bound  total   expected 
cost  with  I  the  inspection  cost  for  each  inspection.     As  it  can  be  seen 
the  rate  of  increase  of  L**  with  the  increase  of  I   is  higher  for  lower 
values  of  I   for  both  cases,   i.e.,  when  the  information  about  two  points  of 
distribution  are  known  or  the  information  about  only  one  point  is  known. 
Table  6  gives   the  values  of  L**  and  I  which  are  plotted  in  Fig.    10  where 
Li*  and  Li*  stand  for  the  values  of  L**  when  two  points  and  one  point  of 
the  distribution  are  known  respectively.     The  third  column  of  Table  6  shows 
the  percent  improvement  in  L**  values  when  additional    point  of  distribution 
is  known  for  different  values  of  I. 

4.7  VARIATION  OF  L**  WITH  a 

Figure  11  shows  the  variation  of  L**  with  a_  the  cost  per  unit  time  of 
undetected  failure.     It  can  be  seen  that  the  rate  of  increase  of  L**  with 
the  increase  of  a_  is  higher  for  lower  values  of  a_  for  both  cases,  i.e., 
when  the  information  about  two  points  of  the  distribution  are  known  or  the 
information  about  only  one  point  is  known.     Table  7  gives   the  values  of  L** 

of  L**  when  two  points  and  one  point  of  the  distribution  are  known  respective 
The  third  column  of  Table   7  shows  the  percent  improvement  in  I**  values  when 
additional   point  of  distribution  is  known  for  different  values  of  a.     The 
value  of  L**  at  a  =  0  is  equal    to  0.8  or  the  cost  of  one  inspection  which 


i  =  20.   $/Un1t  Tiitif 
:  =  8.   $/Inspectior 


0.2         0.3         0.4        0 
x. ,  The  Time  where 
p       =  F(x, )   is  Known 
1   in  Units  of  Time 


x,  For  x„  =  0.6  where 
py  -  F(xJ  -  0.700. 


o,  For  x2  ■  0.5  where 
pv     -  F(xJ  -  0.535. 


Fig.    9.     Variation  of  L**  with  the  Times  x,   and  x2  at  which 
p       and  p       are  Known. 
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Fig.  10.  Variation  of  L**  with  Cost  of  an  Inspection  I. 


Table  6.  Percent  Improvement  in  L**  for  Various  Values  of  I 
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Table  7.     Percent  Improvement  1 


*  for  Various  Values  of  a 


at  least  is  needed  to  be  performed  in  order  to  detect  the  failure  of  the 
system. 

4.8  EXAMPLE  PROBLEM 

Problem:     The  management  of  an  organization  wants  to  introduce  a  new 
deteriorating  item  into  the  system.     There  are  three  deteriorating  items 
that  have  the  same  function  but  have  different  failure  and  cost  character- 
istics.    Because  of  the  shortage  of  investing  capital   the  company  wants  to 
choose  the  alternative  which  has  the  lowest  possible  expected  cost.     The 
characteristics  of  the  items  A,  B  and  C  are  as  follows: 

Item  A  Item  B  Item  C 

1.  Average  inventory  cost,  $/day             140  200  100 

2.  Cost  per  inspection,   $                            56  20  80 

3.  Probability  of  failure  p                  0.180  0.216  0.535 
before  time  x. 

4-     x,   is  equal    to  4th  day  3rd  day  10th  day 

5.     Maximum  life  time  14  days  3  days  16  days. 

Besides  these  the  inspection  can  only  be  performed  at  certain  time  and 
only  once  a  day  for  example  at  3  p.m.  every  day.  The  replenishment  of  the 
item  will  occur  at  the  end  of  periods  equal  to  the  maximum  life  time  of  the 
item.     Which  alternative  should  be  chosen  with  the  above  limited  information? 

Solution:     Since  the  failure  characteristics  of  the  system  are  not 
completely  known  then  the  best  estimate  of  the  costs  would  be  the  maximum 
possible  expected  cost  L**.     In  many  cases   the  systems  show  an  IFR  failure 
characteristic  so  it  is  assumed  that  failure  characteristic  of  the  items  are 
IFR.     Average  inventory  cost  becomes  equivalent  to  the  cost  of  undetected 


failure  a  since  a  failed  item  will   remain  in  the  stock  and  assumes  cost 
until  its  failure  is  detected  through  inspection.     We  assume  that  the  in- 
spection is  perfect  and  it  does  not  take  time  and  does  not  degrade  the 
item.     It  can  be  shown  from  Eq.    (2.4.2)  that  if  L  is  multiplied  by  a  seal 
factor  u,  and  t,   the  time,  multiplied  by  another  scale  factor  v  then  new 
L  called  L'   ■  (u)L  is  the  expected  cost  of  the  system  provided  that 

I1  =  (u)I, 


where  I1  and  a1  are  new  scaled  inspection  and  undetected  fail 
t'max  and  'max  are  the  scaled  and  not  scaled  values  of  the  ma 
time  of  the  system. 

For  Item  A  we  have: 

I'  ■  (40)(1.4)  ■  56 
a'   =  (§)    (35)  =  140 
x'   =  00)04)  ■  4 
'max  =  (10)    (1'4)  =  14- 
So  for  u  ■  40  and  v  =  10,  we  have   I  =  1 .4,  a  =  35. ,  x  = 
and  tmax  =1.4.     From  Table  3  the  value  of  L**  for  above  valu 
parameters   is  L**  =  10.10.     So 


L**1    =  (40)(10.10)   =  404  $/cycle  or  404/14  =>  28.86  $/day. 
For  Item  B  we  have: 

V  «  (100)(0.2)  =  20 
a'   -  (^)(20)  =  200 

x1   =  (10)(0.3)  =  3 
t'max=   (10)(0.8)  •  8. 

So  for  u  =  100  and  v  »  10,  we  have  I  ■  0.2,   a  «  20,  x  =  0.3,   p  «  0.216 
and  tmax  =0.8.     From  Table  6  the  value  of  L**  for  above  values  of  input 
aarameters   is  L**  ■  2.58.     So 

L**'   =  (100)(2.58)  =  258  $/cycle  or  258/8  «  32.25  $/day. 
For  Item  C  we  have 

I'   -  (100)(0.8)  =  80 

a'   =  (Jga)(20)  ■  100 

x'   =  (20)(0.5)  =  10 

t|iax  ■  (20)(0.8)   =  16. 

So  for  u  «  100  and  v  =  20,  we  have  I  =  0.8,   a  =  20,  x  =  0.5,  p  =  0.535 
ind  tmax  »  0.8.     From  Fig.   8  the  value  of  i**  for  above  values  of  input 
larameters   is  L**  «  4.45.     So 


L**' =  (100)(4.45)  •  445  $/cycle  or  445/16  =  27.81  $/day. 
Comparison  between  the  cost  per  day  for  the  three  items  shows 


CONCLUSION 


the  problem  of  partial 
is.  In  this  work  this 

•ogramming  methodology  to 
n'ntenance  cost  of  a 

ng  failure  rate  distribution. 


Minimax  policy  has  been  devised  to  cope 
knowledge  about  the  failure  distribution  of  sys- 
policy  has  been  adopted  together  with  dyn. 
find  the  upper  bound  on  the  optimal  total  expec- 
system  subject  to  deterioration  with  ai 
The  major  difference  with  the  previous 
procedure  and  a  computer  program  has  been  devised  which  finds  the  upper 
bound  cost  not  only  when  one  point  (x,  p)  of  the  failure  distribution  is 
known  but  also  when  two  points  (x, ,  P1 )  and  (x2>  p2)  of  the  failure  distri- 
bution are  known. 

The  basic  findings  can  be  summarized  as: 

1 .  The  convergence  and  accuracy  of  the  upperbound  on  the  expected  optimal 
total  maintenance  cost  depends  on  the  number  of  increments  into  which 
the  feasible  range  of  the  state  variables  c  and  p  are  divided  and  for 
reasonable  amount  of  computer  time  a  fairly  accurate  upper  bound  value 
can  be  obtained. 

2.  Additional  knowledge  about  a  second  point  (x2,  P2)  improves  (decreases) 
the  upper  bound  in  general.  But  it  improves  the  upper  bound  cost 
especially  if  the  second  point  is  closer  to  the  end  point  (T   ,1) 
where  T    is  the  maximum  life  time  of  the  system. 

3.  The  closer  the  two  known  points  (x1 ,  p1 )  and  (x2,  p2)  to  the  end  point 
(T   »  1)  the  more  improved  upper  bound  cost  can  be  obtained. 


n' ch  have  been  found  for  several 
t  all  were  lower  than  upper  bound 


4.  The  upper  bound  cost  seems  to  increase  faster  with  the  increase  in  the 
cost  per  inspection,  I,  at  lower  values  of  I. 

5.  The  upper  bound  cost  seems  to  increase  faster  with  the  increase  in  the 
cost  per  unit  time  of  undetected  failure,  &,   at  lower  values  of  a_. 

6.  Several  optimal  expected  total  costs 
different  IFR  distributions  showed  t 
cost  as  it  could  have  been  expected. 

The  basic  disadvantage  of  the  present  computational  procedure  is  due  to 
the  existance  of  three  state  variables.  In  this  work  the  state  variable, 
time,  has  been  assumed  to  be  discrete  in  order  to  increase  the  computational 
feasibility.  This  assumption  although  limits  the  scope  of  the  applicability 
nonetheless  it  can  be  a  practical  assumption  since  in 
nd  inspection  actions  can  only  be  done  at  dis- 
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works  in  this  area  are  worthwhile  to  be 
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MATHEMATICAL  TRANSFORMATION 
OF  THE  EXPECTED  TOTAL  MAINTENANCE  COST 

As  stated  in  (2.4),  Eq.  (2.4.1),  i.e., 

n-2  Vl 

L(tQ  -  0,  t.,  t, t  -  T)  -  I    J        [I(k+l)+a(t   -x)]f(x)dx 

k=0  t.  K  ' 

t  »T 
+  /  n  [I(n)  +  a(T-x)]f(x)  dx 
tn-l 

n  =  1,2,  ...,  N,       (1) 
can  be  transformed  into  Eq.  (2.4.2),  i.e., 

n-1        rn-l  ln       •, 

I'll     F*(tk)  +  a  I     f*(tkHtM~\)   -  !       F*(t)  dt 

k=0    K    W=0  <   k+i  k   Q        j 

n  =  1,2,  ...,  N.        (2) 
PROOF: 

n-2  Vl  n-2  ,     tk+l 

I     j         [I(k+1)  +  a(t k+1-x)]f(x)  dx  ■  I     \l         I(k+l)f(x)  dx+ 
k=0  tk  K  '  k=0  I  tk 

tk+l  tk+l       ! 

a  /    tk+1f(x)dx  -  a  /    x  f(x)dx  , 


/    I(k+l)f(x)dx  -  Kk+l)(F(tk+1)  -  F(tk)), 


■  KO+IMFU,)   -   F(tQ))  +   I(l+1)(F(t2)   -   F(t1)) 

+  I(2+l)(F(t3)   -  F(t2))  +  ...  +  I(n-2+l) 

(F(tn-l]   "  F(tn-2]) 
-IFCtp)  -  IF^)  -  IF(t2)    ...   -   IF(tn_2) 

■  rtn-DFCt,,.,) 

n-2 
'■I     I     F(tk)  +  Kn-DFC^),  (5) 


\I0{  V^xjdx-a^t^CF^J-FCtk)). 

We  know  that  according  to  the  integeration  by  part  formula 
/     vu'du  =   (vu)^  -  /     v'u  du. 

where  v'    and  u'    are  the  derivatives  of  functions  v  and  u.     Assuming 
u'du  =  f(x)dx  and  v  =  x  and  using  Eq.    (7)  we  can  write 

Vi  t.+1      Vi 

/  xf(x)dx  -   (xF(x))  k+1    -  /  FOOdx 


■  t,+1   F(t       )  -  t  F(t  )  -  /  F(x)  dx, 


n-2  Tc+1  n-2 

I    -a/         xf(x)dx  -  -  a     I     (t.+,F(t       )  -  t.F(t  )) 
k=0  t.  k=0       k+1       k+1  k       k 


Now  substituting  the  equivalent  values  from  Eqs.   (5),   (6)  and  (9) 
in  Eq.    (3)  we  have 

n-2     Vl  n-2 

I     /  [I(k+1)  +  a(t.+1-x)]f(x)dx  =  -   I     I     F(t.)  + 

k=0  t.  k=0 

I(n-1)F(tn_.,)  + 

+  a  /  n"     F(x)dx 
0 

n-2 
-  -   I     [     F(t.)  +  I(n-l)F(t     .) 
k=0         K  n~' 

+  a  f-  ni  F(tk)(t    -  tk)  ♦ 

I     k=0  K   ' 

t     , 
/    "    F(x)dx). 

0  J 


The  second  portion  of  the  right  hand  side  of  Eq.    (1)  can  be  written  in 
the  same  way  as 
t  -T 
/  "     [I(n)  +  a(T-x)]f(x)dx  =  I(n)(F(T)  -  Ftt^))  + 
Vl 

t  -T 
+  a  T{F(T)  -  F(t       ))  +  a  /  "       F(x)  dx 
Vl 

-  a   (TF(T)   -  t^Ft^J) 
=  I(n)(l  -  Htn,}))  * 

a(-  Ftt^Jd  -  tn-1)  + 
/  F(x)dx],  (11) 

Vi  J 

now  substituting  Eq.    (10)  and  Eq.    (11)   in  Eq.    (1)  and  naming  t  =  x  we  have 
L  •  -  I  "]     F(tk)  +  Kn-IJFU^)  +  a[-  "j     F(tk)(tk+1   -  y 

t  _. 

+  /  ""     F(t)dt  +  I(n)(l-F(t     .))  +  a(-F(t     ,)(T  -  t     .)  + 
0  n   i  t        n  i  n   i 

tn=T 
/  "       F(t)dt 

Vi 

*  I  "I     (1   -  F(tk))  +  a(-  ^     F(tk)(tk+1    -  tk)   +  /  "       F(t)dt], 


denoting  F*(t)  as  F*(t)  =  l-F(t)  and  substituting  it  instead  of  1   -  F(t) 
in  Eq.    (12)  we  have 

n-1  f     n-1  VT  , 

L  -  I     I       F*(tJ  +  a  -     I     (1   -  F*tO)(t.+1   -  t.  )  +  /  (1   -  F*(t))dt 

k=0  l    k»0  o  J 


tn«T  t  -T 

►  /         dt  -  /  F*(t)dt] 


I     I      F*(tk)+afj     F*(t.  )(t,+rt.  )  -  /  "       F*(t)dt]. 
k=0  K  >-k=0  K       K+l     k         Q  .  j 
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OPTIMAL    EXPEL 


CALCULATION    CF    THE    JPP 

MAINTENANCE    COST    OF    A    SYSTEM    AlTri    PARTIALLY    K.NC 

DISTRIBUTION 

NUTATIONS 

NB=1 —  CALCULATES    OPPCR    3UUNU    COST    .«HEil    ONLY    ONE 
DISTRIBUTION    iS    k(,u»,\- MODEL    A 

NB=2— -CALCULATES    uPP^R    oCuf-D    LOST    *H<.N    T.,0    PUR 
DISTRIBUTION    ARE    KNOWN — MODEL     B 

PI THE    KNOWN    CUMULATIVE    FAILURE    PROBABILITY     i 

P2 THE    KNOhN    CUMULATIVE    FAILURE    PROBABILITY    ; 

Ai INSPECTION    LOST    PER     INSPECTION 

AA COST    CF    UNDETECTED    FAILURE    PEn    UNIT    TIME 

NP NUMBER    jF     INCREMENTS    INTC    wn  IC.-i    TrtE    FEASIBLE 

FAILURE    PROBABILITY     P    IS    uiVIDEC    MINUS    CMS 

LA NUMBER    OF    R.CREHENTS    INTC    rtHlCh    THE    FEASIBI 

PARAMETER    QF    TRaNFORMED    EXPONENTIAL    OiSTRlBuTIOf 

DIVIDED    MINUS    uNE 

El MAXIMUM    POjSIBLl    NUMBER    OF    I  i\ 


KIT  , 


ujll, 


-Thi 


Xl.iL.' 


JP    TO     TIME    T(I)     *ITH    ?<u),CU)     AND    AT    STai 

RESPECTIVELY 

THIS    PROGRAM    wAS    WRITTEN    BY    K.CU  RuSH-r,  ARAMi 

INDUSTRIAL    ENGINEERING    kA.NSAS    oTATc    DNWE; 

KANSAS,  MAY,  1979 


M-l  ; 


INTEGER  0 1,  02,  £  1,  GO  1,0  DO  l,0D2,ODD2,EE  I 
DIMENSION  RKK( 10,20, ID)  ,  CU0,2u,lG)  ,.1FlU,i3,  1J 
RSB(10,20,10),MK(10,20,lu),P(lu,2C),PK2(10),JJR 


IFlNB-lUCG,o2 
REAG1 5,322)P1 
FORMATliFlO.oi 

02  =  01 

1,3  20 

A  I  ,AA,.\P,c 

415) 

P2  =  P1 

GC    TO    823 

FORMAT <4F 10.6 

PAC=i.E-25 

315) 

REA01  5.600JT1  L 

11 

CIC1--ALOGC1.- 

PU/TILU) 

STEP°llfUR°BOT 

H    MODELS    A 

1^01=01 
002=J2- 
ujQZ=uZ 


82S    PP2=l.-tl.-PJ 
IF(PP2>3,4,4 

3  DP=PP1/P,M1 

Gu    TO    6 

4  DP=(PP1-PP2), 


go   ro 

9    P(L2 

1C     lui'-ll 


3=2, M?l 

3)=?(L2,L3-1)-0P 

1)100,3,386 


-ll.-Pl)*b; 
-EXP(-CiCX» 


P1-?P2)/PN1 
l)  =  i.-PAC 

Lr  =  2,NPi 


13  P1L 
P(L 

14  CONTINUE 

STEP    2    F-CR    BC 

OC=ClCl/ALi 
C(l, 1,11=0. 

2  7   C(  1,  l,iO  =  C(  1, 


16   L 


1(1 


L10 


17   COi 

IFINS-1U 

33  00=10  102- 
t 101 ,1  ,1) 
00    23    r.  =  2 

28    CID1,  l,KJ: 

C (01,1 ,La 
00    20    LS*I 


00=(C2-C 
OC    19    K  = 


2  1    0  I  02,  1,K)  = 
0(02,  1,LA 
00    It    LI 
00    512    Ll. 
IF  (LI  2-1) 

707  Cl—ALOGU 
GO    TO     70  i 

708  Cl  =  -ALOGl 

709  ClLli ,L12 


;,L10)I  )/(T(D2)-T(L' 


=0li_9,H0,is-l)+00 


02* i  n  ooo 2) -r  (0  2)  i 


i  T(  Ou02i-r(  02) 


=  1,NP1 
00,7C7,70ci 
AC/U.-P2)  )/  IT 

l.-P(LU,L12)), 
11  =  01 
512,24,24 


7  00  CtL13,L14,K)*-AL0U((l.-PPl-Ut 

113+11-TIL13)) 
610    C(Llj,Llt,LAi=-ALCo-(PAO/(  l.-Pi 

STEPS    3i-,    FCR    BOTH    ,100cLS    AiO 

I  =  D1 

RKK101, l,2)=Ai+AA*T(01)-AA*?([ 

00    25    \=2,LA 

SsKIOl,  l,KI»ftKK(Dl,  1,2) 


/ITILUtU-TIL 


R*K(02,1,K 

MF102,] 
MT(D2,J 
MK<02,1  ,K)=1 
00    Z<i     1=0002, tl 


2S   CONTINUE 
570    WRITE(6,8 

o04   FCRMATUX 


,     dUTH     .-'.OCElS     Ai3     r 


HRITE(6,32)H 

32  FCRMAT(2X,'M=' 
00  30  1=2, 0C1 
00    30    J»2,NP1 


3C   CONTINUE 
STEP    6    FOf 
loo    DC    31    II=f 


IL=1 

[  [»H 

42    IF  UJh 


145    -IT 

Ut    I  I' 


2(H)  --.J.-< 
Fil i  1  )  =  J 
<  =  J*  +  1 
[JR-NPU 
R( II  )  =  JR 

3=PK2 (11 

2=11 

=  11  +  1 

iPK2<m 

=  II  +  L 

v3=PK2llI 

(Pl<2(  11  )- 

2  =  11 

=  11  +  1 

I  11-001) 
[RKK1D1, 
KID1,1,K 
(Ol.l.K) 

IF(Na-l)  100  ,660,680 

180  M.M=0J2 
JVR=2 
0:5    TU    861 


5t   PPK2=dOOQ. 
JR=2 
L=l 

FP=C(  I,  j,  1J*(  TIul)-T<  II)  ) 
IF(FP-10:).)<io3,369,3o9 
069    ?PP1  =  0.0 
GO    TO    49 
368    PPPl*i.-ll.-P(01,l))*exP(C< 
49    IFIPPP1-PI  I  I,  JK)  )  50,50,51 
5  0    IFIKKKI  II,jK,LA)-4D30.)3ol,5J,53 
3  81    PPK1=(  U+AA*(  III  )-I  i  I  1)  )  )*(  1.-PII  £,.„ 
lCiII,JR,UAJ+(P<Dl,ll-P(I,J))/CtIf  J,l 
IFIL-1  )  100,68  ,  52 
od    L  =  2 

PPK2»PPK1 


T(01)-T(  II  J  J  ) 


IFUR-NPU49, 49,51 
52    IF(PPK.2-PPK11  34,d3,53 
5  3    Jk=JR+i 


IF(JK-NP1>4 

9,49 

.51 

51    IFUL-IHOO 

i55  , 

56 

PPK3=PPK2 

61   MFM2=MFM1 

0.J6 

MTM1»U 

62    11=11+1 

I  Ft [I-UOXJS 

,57 

53    IFIPPK2-PPK 
uC    I  1  =  11  +  1 

3)59 

,60,60 

IFIH-OD115 

59   PPK3  =  PPK2 

6,56 

,57 

IFIPPK2-700 

0.  16 

64    MFM2=MFM1 

.-1T.-1 1  =  1  I 

65    11=11+1 

IFUI-OC115 

,57 

57    IFIRKKtl ,J, 

li-J 

PK3)b7Q 

67   RKKi  I  ,J,1)= 

PPK3 

rtF(I,J,l)=M 

F;12 

7  2    IF  (.-UK  I  II  ,Jk,LA)-40J0.)332,7-,,7  4 

6  2    PPK1=(AI+AA*ITIC2)-T(  I  I  )  )  )  *  U  .-?  I  1 1  , 

ILIII  ,  jR,LA)H?lul,l)-P(02,l)  1/CIC2H 
IFlL-1)  100,  70,  71 

7  0    L  =  2 


1)72,72, 75 


7  3    ?P*2=P 

JS=jK+ 

IF  (JR-I 

75    IFiLL- 


Si 

»*K3=PPK2 

37 

|FM2»rtFMl 

ac 
a6 

11=11+1 

IFUI-OOU  79,  79,  30 

IF<RKK{02,1,2)-PPK3)85,8 

HKM02,lr2J«PPK3 

*F(D2,l,2)«MFrt2 

Mnu2,i  ,zi~mm 

1K(D2ti,21-« 

jq  37  k=2,la 

37 

1FI02, 1,K»*MF»02, i,2) 
0102,1  ,K)»MT(02»1.21 

*U(J2,  1 .  K  )  =  ,1 

35 

DU    135    I  =  OC02,E1 
JC    105    J»l ,NPl 
I  I=M 

99 

lL=1 
?PK2=6000. 

/CIC2+{P(02, 


92,92,95 
>K1)93,S4,< 

92,92,-j5 


-7900.1390,391 


I F ( I I-UL 
9  7  tF(PPiO- 
102    i  1  =  11+1 


9  2   MFM2«r 

HTMX'Il 
93    I  1  =  1 1  - 


535    IFU-E1  )i05  ,530,100 

105  CONTINUE 

53  0   KRITE(6,605JPPKJ,P(MTM1,MFM2 

3  05   FbRHAT(4X, 'R0='  ,H2.6,2X,'  J  = 

STEP    7    FOR    BOTH    MCOELS    ASo 

lFlH-OOU180,a45,  100 

645    IFlNd-1 )100, 353,21a 

STtP     5    FOR    dCTH    MODELS    AtS    M 
laO   M=M+1 

*RiTE(6,  125M 
125   FCRMAT(2X,'M='  ,15) 

00    129    t«M,OUl 

00    129    0  =  2,  ,V1 

Il  =  iM-l 

12=1-1 

00    111     11=11,  12 

PK2U  I)  =6000. 

PPPl*i.-(l.-P(i,Jl)*EXP(C(  I, 

00    109    0R  =  2,1mP1 

•PtI  IiJRJJ10S.ua, 110 


v-ONTI 


GO  TO  111 
110  JJR(U)  =  JR 
1  1  1  CONTINUE 

00  129    ,<=2,LA 
IF(I-0Ql)52j,524,iC0 

524    IFU-LAJ  129,523,100 

1  1=M-1 
12=1-1 

JC113  i  1  =  1  1  ,  12 

IF (JJRI II i-NP) 112 ,54i 
112  PPP2=l.-(l.-pj  [,J))*E 

J*-=JJR(I  I  ) 
7  IF(PPP2-P(1I,JRJ  )  114, 
4  C1=-AL0G((1.-P(  I,  J)  u 

IF(M-3)10Q,171,172 
2  00  173  RP=2,LA 


71    lFlhUill  ,JR,L*)-40J0.) 
S5    ?PKl=(AI*riA*lT(i)-T(  II) 

l  +  RKKdl  ,JR,LA) 
21    1FIPK2I  IIJ- 7000.  )  liO.li: 


ISO    IF(PK.2(I1)-PP,U)U: 


men  m=jR 

JR=JR+1 

112 

IM  jr<-NPl)  117,  117,113 

jjrii  i)»jr 

547 

IF(LL-1  1100,119,120 

PPK3=PK2(Il) 

IF(PK2lli)-700J.H21,lio,Ua 

121 

HFM2=MMF1UI) 

GO    TO    113 

120 

IFIPPK3-PK2UI)  )  118,113,  123 

123 

PPK3=FK2(I1  ) 

IF1PK2I  IH-7O0J.)  124, 113, lid 

<FV=.'.;<KI  1  i  ! 

MTMl«II 

CONTINUE 

RS8(I,J,K)=PPK3 

1FIPPK3-7000.  )  127, 120, Ud 

127 

«-<ITt(6,1261Rd6(I,J,K),PHT,ll, 

126 

F0R,MAT12X,'R='  ,F12.6,2X,  '  F- >     F 

,F12.c,2X,'C=' ,Fl4.ol 

126 

"'1  Tc  (6,13u)KJ,a(  I  ,J,k)  ,  I  ,J 

ill 

F0RMAT(2X,'R=.»,fl2.6,2x,.l»i,i 

JO    730    1=M,DJ1 

JO    73  0    J=2,NP1 

tFU-uQU  731,732,100 

73  2 

IF  IK-LA  1730, 731, 100 

DG  ZZO  *  =  2  ,LA 
1VI-MFID1,  l.KI 
MV2«MT(01,l,K) 

1P(  I,  J)  ,CI!,J,<|  'U  '"'    ' 

2.1    FCRMAT(2X,'R»"fH2.o,2x,'l 

l',I  =  ,.I5,2X,-J=.,F12.6,2X, 


DO    223    J=2,,V1 
RKKUfJ,LA)»RKKU,J 
MV1=MF<  I  ,„,l) 
1V2»MT(I,J,1J 
WRITe<6,224)KKK<  1  ,J, 


STEP    12    FOR   MODEL    B 
:   KRUE(6»225)M 


UU    226    1=11,002 
OJ    22  &    J*2,NP1 
I3»M+0Ui-l 
12=1-1 

JO    227    11=13,12 
PK2U  11*8000. 
IFCII-01)100,227,313 
313   PPPl»l.-U.-PU,J))*EXP<C(l,J; 
JO    229    JR=2,NP1 

229  CONTINUE 
JJR(I I J  = , 4 P 

Gu  TO  227 
210  jjrii  I  I=jr 
227   CONTINUE 

JO    226    K  =  2,LA 
IF(I-DD2)550, 551, 1 J a 

551    IF  (,\-LA)226  ,550,100 
55C    LL*1 


232   PPP2=1.-<1.-PU,J). 

JR»JJRUl) 
249    lFiPPP2-P<U,JR})2 
236    Ci=-ALOGI (  l.-H  1, Ji 

IFIM-2U00, 237,236 


1239,240,240 
H300.)H02,248,24a 


230   JO    Z^Z    K.?=2,LA 
IFICIOl.l.KP)-! 
2*2   CONTINUE 


1.1)-P(I,J1  l/C(  i,J 


■NKi  = 


(II 


IF(PK2UI)-7QOt; 

254  MFM2*«Mfi(  1  I) 
MTM1  =  II 

GO    TO    223 
253    IF(PKK3-PK2(U))228,228,255 

255  PKK3=PK2UI ) 
IF(PK2(  I  U-7000. 1256,228,223 


H TM 1=  II 


I  I) 


00  TC    226 

2  =  S  »fRITc(b,260)R8a(  I  ,J, 
260  FQKMATI2X, • R= ' ,fi2.6 
22t   CONTINUE 

1  l«DDl+M 

JCj    735    1  =  11,002 

00    735    J-2.NP1 

OG    735    K=2,LA 

IF([-Q02;73o,737, 133 
737    IFIK-LAJ  735,736,100 
736    RKK(I,JtKJ«R88(I,J,K; 
735  COMTINUE 


STEPS    10,; 


■    FCR    , 


U«M+D01 

30    301    Ii  =  I  1,032 

IF(II-Oi)100,301,303 


30    305    1 1  =  I  1,002 
IF(II-01U0O,3O6,3C7 

I    lF(JJRim-NPlJ308,JOd,309 


07    PPK1  = (A 
ICdl,  JR 

GO    TO    3 


•PlIl.JK)  1310,3 
I  ,JK,LA)-<»000. 

+  AA»IT[D2)-Il  I 


*(  l.-P(Ql,  1) 


.1321.321.311 


-7000.  132*.  3i.i,  JJ3 


23    IFIPK2UI] 


326   ,1TM1  =  II 

MF,M2=.MMF1(  I  I 
305   CONTINUE 

IMKImsI  02,1,1 


PP<3)325,3u5,305 
7D00.)326, 305,305 


DO    330    0=1 

,NP1 

U*M+0Q1 

DO    331     11= 

V0J2 

PPK2»8OO0. 

00,332,^33 

Jrl  =  2 

FP  =  C(  I.J.L 
IFlFP-iOO. 

cGJ.eOS.oOV 

ppp>i!-u 

-P(  J2.1)  )*EXPICI 

IFUKKUI, 

J  A,  LAI— +000.  K'.ii  , 

PPK1»(AI*A 
tlll.JR.LA 

4*(T< I 1-TIII 1 1 1*1 
♦  I P< 02,1  I -PI  I, J) 

IF(L-XU00,336,337 

PPK2=PPK1 

RFH1=JR 

JR=JR+1 

If  (JR-NPiJ333t33ai335 

337    IF(PPi<2-PPKi)340,341, 1' 

IF(JR-NPI)338,338,335 


PPK3=PPK2 

I(M?PK2-7JC0.J345,33i,331 
3  45    MFM2=MFM1 


:.J,l)-HPK3)532,532,o03 


603  RKK! I ,J,1)=PPK3 
■'•II- (I  ,J,1)=,«F.12 
MT(I,J,1»«MTM1 
MK(I ,J,1)=M 

532    IF(I-EU330,3i2,100 


IF  (,1-02*01)350,350  ,353 

STEPS    aS9    PCR    MUoEL    A,     STEI 


362   FGR.4AT1 


(  I.  J,l  i,  I  ,t>  II  , 


=  ■  ,t-12.6,2X,,i 
*30il00 


17    FGS    MOOEL 


3  71  UG  414  Jk=1,: 
IFIRKMIl.JR, 

tlC  PKIU»IAH-AA*( 
IF(PKK2-7u00. 

+  13    PKK2=PKM 


416    PnK.3  =  PKK2 


LoJ,<fi  = 


-/GOO.  )4lb,372,372 
£    •IFM2=:1F;U 

MTH1=H 

GC    TG    J72 
7    IFIPKK3-PKK2J372.372.419 
S    PKK3  =  ?.\K2 

IF(P,NK2-roOC.)42C,372,372 
3   MFM2=MFM1 

,'-1TMl=lI 
2    CONTINUE 

i-.KITt  I6,300)P<.<3  ,  PUT,-tl,,4Ff 
:    FGRtfAT(4X,  'RC  =  '  ,F12.o,2X,'. 

IF(PKK3-Ri 


423    CONTINUE 


IFtM-EH-021431,' 
STtP     12    FCk    MOOI 
3  1    WRITE  16, 470  ),«l 


12«M+0D2-1 

13=1-1 

UO    435    11  =  1, 

i.  13 

?KK2=d0J0. 

IF  U  1-02)10) 

►36    LL  =  2 

JU    SOI    K»2,l 

72 

,•.3  7 

>01  CONTINUE 

503    PKK3=(AI+AA' 

i-tic; 

1C(  I  ,J,1)+K\i 

LI  0 2 ,  ] 

*FM2*1 

MTM1=02 

Gb    TC    435 

f37    PPPl=l.-ll.- 

-P  (  I ,  ] 

.  ))*EJ 

IF(JR-NPl) 441,441,459 

-3S    JJV=Jrt 

-50    IF <M-2)100,690,620 

QL    o2  1    ,<P=2,LA 
IFIC(U,JJV,KP)-C  1)621, 622, 622 

o21   CCNTINUt 

622    IFUKiCUI  ,  J.v,<P] -4000. )a23, -45  ,445 
6  23    PXK1=(AI+-AA*(T(  1  )-T(  11)  J  )*[  l.-PUI,J. 
1/C1+RKMI1,JJV,KP) 
Go    TO     625 
690    I FIRKK ( II , J JV, LA  1-4000. 1443,445, 445 
443  CI— AU)GtPAC/(l.-PUI  ,JJV11  I/ITUI-T) 
PKK1=(AI+AA*(T(I)-T<II)))*(1.-PU1,J. 
1)  l/ol+kKiUII.JJV.lA) 
o25    IFIPKK2-7000.  1  447  ,44b  ,44d 


I  1-  (JjV-i\Pl  1450,45  J,  -,59 
4,7    IF  (Pim^-PmU  1-52,456,4: 

-:>2    Pnn2=?i<K1 


IMPKK.3-7000, 


I    lr(PKK.3-PKK2)4j5,435,4fc-. 
(    Pr,N3  =  PK.K2 

[F{PKK3- 7000.) 465,435,435 


,J,LA)=PKK3 
K3-7000.)  til, 473, 4: 
)RKK(  I, J, LA) 

R««  ,F12.6,2,< 


473    ml 

475   FC 

705    JLi 


710  PPP2=0. 
3G    TO     711 

712    PPP2=  l.-(  1 

711  JR-JJRUl) 

PPP2-P1. 

M-21100 
CI  II 


6fc3, 663,64b 

0,667,<L6o 


PKK3 

10.  )b72,673,o73 
>)R3Bi  I  ,J,M  ,PC 


DO  740  1= 
DO  740  J» 
00  740  K= 
IFU-ttl) 


'40  CCNTiNUE 

STEP  13  FOr.  MODEL 

30  TU  430 

STEP  14  FOR  MOOtL 

ouUAJ  E^PEuTEO  Tl' 


APPENDIX  C 

COMPUTER  PROGRAM  FOR  CALCULATING  THE  OPTIMAL  EXPECTED  TOTAL 
COST  WITH  COMPLETE  INFORMATION  ABOUT  FAILURE  DISTRIBUTION  OF 
THE  SYSTEM. 


calculatil: 


MOTATIQNS 

AI INSPECTION    ^bST    PEri    INSPECTION 

Art  —  CuST    ur    u.'.jCTelTcJ    FAlLJRt    P  br.    JN 

1(TF) :1AXHIM    Lil-E    TIME    CF    THE    S^SIl 

KK(I)  ,K3(I)  —  -oPTHJri  EXPECTED  rtAINTcN. 
1TIME  T(I)  ANO  AT  STAGt  M-l  A.MO  f  RE5PH 
THIS  PROGRAM  WAS  »RITT£>;  3Y  <OUROSH-KM 
.ANSAS    STATE 


7C    Tne 
lEPARTMENT 


iMSAS, 


979 


1  FORMAT12F10-* 

INPUT     OATA     Ac 
DO    2     1=1,1 

*0   FORrfATOfi 

2  CONTINUE 


(I  )  ,C( I)  ,T( I) 


CaLCJLATES    R;UTM     FCR    tf>2 

3C  uu    13    1 l«M,TFl 

UC    31    LL=  I  I  ,TF1 
31    V=V+(P(LL)-P(LL-H))/C(LL  +  1) 

PK»(AI*AA*<  r(TF)-T(IIl))*U.-P<II  ))  + 

COMPARES    inMTFI     FOR    CLCcR    ANJ    CURRcN 
THE    3,-IALlER    VALUE    ANJ    QioCARJi    TriE    L 

14   SMTFI  =  P< 
MT*II 

1-3   CONTINUE 


19  CONTINUE 

»RITE(6,2 

id  CCNTINUE 

'   *  'Pl"-a 

00  23  LC= 

2  3  RiULCi  =  Ru 

GO  TO  30 

PRINTS  Tn 

E  OPTIMAL  E 

NUM8ER  DP 

INSPECTION 

vC  *RITEt&,9 
;Y    FORMATUX 

1)RMTF),MT 

cO     rOTAL    COST    ,<MI)     A.-lO    OPTIMAL 
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A  considerable  work  has  been  done  in  the  field  of  maintenance,  but 
except  for  a  few  almost  all  of  this  work  has  been  based  on  the  complete 
knowledge  about  the  characteristics  of  the  deteriorating  system  or  equipment 
especially  about  the  failure  characteristics  of  it.     But  in  real  situations 
the  complete  information  about  failure  distribution  of  a  system  is  rarely 
available  especially  when  a  new  one  is  introduced  into  the  existing  system. 
Minimax  policy  has  been  devised  already  to  cope  with  this  problem.     Basically 
it  gives  the  best  maintenance  or  inspection  policy  or  timing  when  the  in- 
formation about  the  system  failure  distribution  is  incomplete.     It  utilizes 
the  information  which  is  available,  to  minimize  the  total   cost  and  at  the 
same  time  maximizes  the  total  cost  with  respect  to  all   possible  and  feasible 
values  of  the  unknown  portion  of  the  failure  characteristics.     The  present 
work  utilizes  the  information  about  one  and  two  points  of  increasing  failure 
rate  distributions  to  find  an  upper  bound  for  the  optimal  expected  total 
maintenance  cost.     The  basic  difference  of  this  work  with  previous  work  in 
this  area  is  that  a  procedure  and  computer  program  has  been  divised  which 
utilizes  the  available  information  not  only  about  one  point  of  the  failure 
distribution  but  also  searches  for  improved  upper  bound  for  the  total   cost 
when  information  about  two  points  is  available. 

The  comparison  has  been  made  on  the  basis  of  the  available  knowledge 
about  a  failure  distribution  between  the  upper  bound  total   costs  and  the 
value  of  information  has  been  discussed  through  an  example.     The  yariation 
of  the  upper  bound  cost  with  changes  in  the  cost  per  inspection  and  the 
cost  per  unit  time  of  undetected  failure  has  been  discussed  through  an 
example.      It  has  been  found  basically  that  having  information  about  two 


points  of  the  distribution  improves  (decreases)  the  upper  bound  total  cost. 
Also  it  has  been  found  that  the  closer  the  known  time  of  failure  probability 
to  the  maximum  life  time,  the  higher  the  value  of  information  and  the  lower 
the  upper  bound  would  be.  It  has  been  found  that  the  upperbound  total  cost 
is  more  sensitive  to  changes  of  the  inspection  and  undetected  failure  costs 
at  lower  values  of  these  costs.  Computer  program  has  been  written  and  used 
for  several  different  increasing  failure  rate  distributions  to  find  the 
optimal  total  costs  and  policies.  These  optimal  costs  were  compared  with 
the  upper  bound  total  costs  and  all  were  lower  than  the  upper  bound  cost  as 
it  could  have  been  expected.  Finally  an  example  problem  has  been  worked  out 
to  illustrate  the  application  of  the  results  of  this  work  in  industry.  There 
is  still  a  lot  of  study  needed  in  the  field  of  maintenance  with  partial 
knowledge  to  find  efficient  and  computationally  feasible  methods  to  optimally 
utilize  the  available  partial  knowledge  about  the  system  failure  character- 
istics for  finding  improved  upper  bounds. 


